ANALYTICAL STUDY OF 
INITIALLY BENT STEEL PIPE PILE 
IN LAYERED SOILS 


A Thesis Submitted 

in Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


By 

KAVURI SIVA RAMA PRASAD 


to the 

DEPARTMENT OF CIVIL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

SEPTEMBER, 1981 



C£- PtA'ANR 



- 5 MAY 1982 



CERTIFICATE 


This is to certify that this work "Analytical 
Study of Initially Bent Steel Pipe Pile in Layered 
Soils," by Mr. Kavuri Siva Rama Prasad, has been carried 
out under my supervision and it has not been submitted 
else where for a degree. 



(Umesh Day ah/ 

Assistant Professor 
Dept, of Civil Engineering 
Indian Institute of Technology 
Kanpur - 208016 


September, 1981 



ACKNOWLEDGEMENTS 


I take this opportunity to express my deep sense 
of gratitude to my thesis supervisor, Dr. Umesh Dayal 
for his excellent guidance, tangible criticisms, encourage- 
ment for original thinking * 

I express my special thanks to Dr. Kankanalapalli 
Kurma Rao and Mr. Bavisetty Satyanarayana Murthy for 
helping me in various phases of this thesis. 

I would like to thank the management of Computer 
Centre, for providing me excellent services. 

It will be unfair if I do not express my thanks to 
Dr. O.P. Misra, Health Center, IIT Kanpur without whose 
care of my health this work would not have been completed. 

I extend my whole hearted thanks to one and all 

whose presence made my stay at IIT Kanpur a pleasant one. 

* 

Last, but not the least, I thank Mr. O.P. Sharma 
for his neat tracings and Mr. J.K. Misra for his careful 
and fast typing. 


(Kavuri Siva Rama Prasad) 



iv 


CONTENTS 

Chapter Page 


LIST OF FIGURES vi 

NOTATION AND SYMBOLS viii 

ABSTRACT xi 

I* INTRODUCTION 1 

II. LITERATURE REVIEW 7 

2.1 General 7 

2.2 Field and Laboratory Observations 7 

2.3 Performance of Bent Piles and Design 

Recommendations 12 

2.4 Theoretical Analysis 13 

III. AXIAL LOAD CAPACITY OF INITIALLY BENT END 

BEARING PILES IN LAYERED SOIL 18 

3.1 General 18 

3.2 Capacity of Straight Steel Pipe Pile 19 

3.3 Capacity of Bent Pile 21 

3.4 Statement of the Problem 25 

3.5 Assumptions and the Governing 

Differential Equation 26 

3.6 Closed Form Solution 28 

3.7 Other Possible Shapes 40 

3.8 Initial Curvature and Capacity of 

Bent Piles 44 

3.9 Results and Discussion 45 

IV. INITIALLY BENT FRICTION PILES IN LAYERED 

AND N0N-H0M0GENE0US SOILS 54 

4.1 General 54 

4.2 Finite Difference Technique 55 

4*3 Differential Equation for Initially 

Bent Friction Pile 


57 



V 


Chapter 





Page 


4.4 

Analysis of Bent Friction 
Cohesive Soils 

Piles 

in 

58 


4.5 

Analysis of Bent Friction 
Cohesionless Soils 

Piles 

in 

69 


4.6 

Analysis of Bent Friction 
Non-Horn ogeneous Soils 

Piles 

in 

74 


4.7 

Load Carrying Capacity of 
Friction Piles 

Bent 


76 


4.8 

Analysis of Bent Piles Subjected to 
Down Drag 

77 


4.9 

Results and Discussion 



78 

V. 

SUMMARY AND CONCLUSIONS 



86 


REFERENCES 



93 



vi 


Figure 

1.1 

1.2 

2.1 


3.1 

3.2 

3.3 

3.4 

3.5 

3.6 

3.7 

3.8 

3.9 

3.10 

3.11 

3.12 

3.13 

3.14 

3.15 


ft 


LIST OF FIGURES 


Classification of piles based on their initial 
shape. 

Winkler Model 

Deflection components for 14 BP 89 and 14 BP 79 
piles measured approximately normal to and parallel 
to even K axis of each pile. 

Terminology 

Variation of soil modulus with depth. 

Idealized variation of soil modulus. 

Pile boundary conditions considered, respective 

initial shape, BM and SF diagrams - 

Effect of a on the variation with depth. 

Effect of a on the variation of with depth. 
Effect of C/T on the variation of w ^ with depth. 
Effect of C/T on the variation of with depth. 
Effect of R q on the variation of w^ with depth. 
Effect of R q on the variation of with depth. 
Capacity of bent pile vs. z m . 

Capacity of bent pile vs. C/T. 

Capacity of bent pile vs. C/L. 

Capacity of bent pile vs. T/r. 

Effect of R q and z^ on bent pile capacities for 
different end conditions. 



vii 


Figure 

3.16 

3.17 

3.18 

3.19 

3.20 
4.1 


4.2 

4.3 

4.4 

4.5 

4.6 


4.7 


4.8 


4.9 

4.10 

4.11 

4.12 


4.13 


Percentage increase in deflection vs.C/T. 

Percentage increase in bending moment vs. C/T. 

Typical shapes for different end conditions. 

Capacity of bent pile vs.C/T for different shapes. 

Capacity of bent pile vs.R/T 

Variation of soil modulus and axial load in pile 
with depth. 

Bent pile divided into segments. 

Approximated variation of soil modulus. 

Effect of friction on w ^ and ^2 ■‘■ n c °kesion soils. 

Effect of friction on and i n cohesionless soils. 
Capacity of bent pile vs. z ffi in cohesive soils. 

Capacity of bent pile vs. z m in cohesionless soils. 
Capacity of bent pile vs. C/T in cohesive soils. 

Capacity of bent pile vs. C/T in cohesionless soils. 

Capacity of bent pile vs T/r in cohesive soils. 

Capacity of bent pile vs. R Q , z^ and 9-^ in cohesive soils. 

Capacity of bent pile vs. R Q » z-^ and 9^ in cohesionless 

soil. 

Effect of polynomial variation of soil modulus on 
capacities. 



NOTATIONS AnJ SYMBOLS 


&0 * 9 a 2 3 

a 3» a 4» a 5 


i Constants 


A 

A(I) 

b 0 > 

b 2* b 3 
B(I) 

C 

c,, Cp, 

C 1,C 2 

C u 

d 

D 


Area of Pile 

Recurrence Coefficients at node i 
Constants 

Recurrence Coefficient at node i 
Initial pile offset of bent pile 

Constants 

Undrained shear strength 

Internal diameter of steel pipo -pile 

External diameter of steel pipe pile 


Constants 

Young’ s Modulus of steel. 

Shear Force at top and bottom of pile 

FA(I),FB(I) 

F€(l) Coefficients at node i 

h Increment length 

i Node I 

I Moment of inertia 

k Coefficient of subgrade reaction 

K ■ Soil Modulus or Subgrade modulus 


^ 1 ^ 2 * 

d i*d 2 


E 


F 1» F 2 



ix 


K 


z 


L 

L 1 


M l- M 2 

n 


^ht'^b 


cr 


*1» q 2 


r 


r 


© 


R 


R o 

RHS* (I) 

s 1 ,s 2 ,s 3 


Dimensionless depth at which initial shear 

force is zero 

Length of the pile 

Depth or thickness of top layer 

Constants 

Initial and additional bending moments 

Number" of segments or exponent in the polynomial 

variation of soil modulus 

The coefficients of modulus variation for the top 

and bottom layers 

Load applied on the bent pile 

Critical buckling load 

Load in the pile at any depth y 

Initial and additional soil reactions on the bent 
pile 

Radius of gyration 

Non-dimensional soil modulus a "t the pile head in 
non-homogeneous soil 

Radius of curvature of initial- shape of the 

bent pile 

Layer coefficient 

Right hand side of differential equation at node i 
Initial flexural stress, additional flexural stress 
and axial stress respectively 



X 


T 


u,u 1 # u 2 


u 2t» u 2b 


w, w lf w ^ 


( W 2>i 
w 2t* w 2b 


W2(I) 

y 

z 


z 


1 


m 

z(l) 

a 


a 


c 


<Pl» <P 2 

6 


Relative stiffness of the pile 

Total deflection, initial deflection and additional 
deflection respectively 

Additional deflections in top and bottom layers 
respectively 

Nondimens ional total deflection, initial deflection 
and additional deflection respectively. 
Dimensionless additional deflection at node i 
Dimensionless additional deflections in top and 
bottom layers 

Dimensionless additional deflection at node i 
Depth along the pile 
Non-dimensional depth coefficient 
Non-dimensional thickness of top layer 
Non-dimensional length of the pile 
Depth at node i 

Non-dimensional load on the pile 
Load carrying capacity of straight pile 
Skin friction factor for homogeneous or non- 
homogeneous soils 

Skin friction factors at the top and bottom layers 
Non-dimensional skin friction parameters 
Initial slope at the bottom of Fix-Pin pile 
Initial slopes at the top and bottom ends of 
pin-pin pile 



xi 


ABSTRACT 


The observations of pile bending during driving have 
been reported as early as in 1930. Since then many investi- 
gators have discussed the bearing capacity of bent piles. 
Response of the initially bent pile to the externally 
applied load is one of the complex soil structure inter- 
action problem. Review of ^literature indicates that no 
systematic analytical work has been done on initially bent 
piles in non-homogeneous soil which is commonly encountered 
in practice. An effort is made in this investigation to 
analyze the behaviour of initially bent steel pipe .piles 
in layered soil using Winkler s hypothesis. It is assumed 
that the coefficient of subgrade reaction for each layer 
is constant and the ratio of the subgrade reactions of top 
layer to bottom layer is represented by layer coefficients. 
Polynomials satisfying both geometric and natural boundary 
Conditions are selected to represent the initial shapes of 
the bent piles. Three end conditions, fixed— fixed, fixed- 
pinned, and pinned-pinned are considered in this analysis. 

The load carrying capacities of straight piles are calculated 
according to the AISC specification for the Design, Fabrication 
and Erection of Structural Steel for Buildings and the load 
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carrying capacity of the bent piles are expressed as a 
percentage load carrying capacity of the straight piles. 

The effect of initial shape on the load carrying capacity 
of bent piles is studied assuming number of polynomials 

for each set of end conditions. The effect of initial 

% 

shape on the load carrying capacities of bent piles is 
shown to be most significant. The effect of layering on 
the capacities of bent piles is found to be negligibly 
small. 

The influence of side friction on the capacities of 
initially bent piles is investigated for layered cohesive 
and cohesionless soils using finite difference technique 
and the method suggested by Glesser (1953). Down drag effects 
are investigated for cohesive soil. The effect of polyno- 
mial variation of subgrade reaction on the capacities of 
initially bent piles is also investigated. The increase in 
the capacities of bent piles due to soil friction is shown 
to be significant in case of cohesive soils especially for 
smaller pile offsets and longer lengths. The down drag effects 
are found to be severe in reducing the capacities of bent 
piles. The polynomial variation of soil modulus has a minor 
effect on the capacities of bent piles. 



CHAPTER I 


INTRODUCTION 


Pile foundations are generally resorted to -when 
either the subsoil is poor, the loads to be carried are heavy 
or the settlement becomes excessive. A Pile transmits the 
load from the structure to the soil either by point resistance 
or side friction or both. The design and construction of 
pile foundation is still an art and a challenge to the designer. 

In designing the pile foundation it is generally assumed 
that the pile member beneath the ground is in perfect condi- 
tion i.e. straight and free from initial bending stress. But 
perfection is rare and the pile foundations are no exception 
to this. Engineers associated with piling over the years 
have cone to recognize the fact that not all the piles are 
driven straight, no matter how carefully it is driven. To 
account this very often driving tolerance such as less than 
1 in 80 or similar recommendation is specified for construction . 
Awareness of the problem of piles not driven straight arose 
with the extensive use of steel piles in Norway in the 1930s. 
Since then, many cases have been reported as listed in 
Chapter II. 



2 


The classification of piles based on the initial shape 
is shown in Fig. 1.1. When a pile is smoothly and uniformly 
curved without any sharp kinks, it is called bent pile. The 
bending can., be in one or more planes. A pile whose axis 
shows some deviation from the perfect straight line before 

» 

loading is known as imperfect pile. A pile may have kink or 
develop a very sharp curvature locally instead of smooth and 
uniform curvature. Such a pile is called doglegged pile. 

In practice, the reasons for piles not driving straight 
may be one or more of the followings (1) lack of straightness 
in the pile, (2) slight ecentricities of load application 
during driving, (3) asymmetric failure pattern in the soil at 
the pile tip, ( 4 ) the longitudinal and transverse vibrations 
imparted to the pile during driving, (5) the presence of 
adjacent piles, other obstructions and near by excavation, 

(6) variations in soil density, (7) slope of bed rock, (8) splic- 
ing of various segments of pile and (9) other considerations 
such as slackness in the guides as they determine the direction 
of entry of the pile into the ground. 

The pile bending phenomena during driving has important 
practical implications if (a) piles are members of a group and 
a pile which has diviated affects the driving of other members 
of the group, (b) The structural integrity of the pile is impa- 
ired due to large deviation during insertion. (c) Pilesdesigned 
to go through upper layers of soft soil and bear on firm 
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underlying strata due to large deviations, (d) Due to large 
deviations* the application of the working load produces 
additional deflections of the pile which are sufficient to 
cause yield of the soil, the carrying capacity of the pile 
might then he impaired. 

The capacity of the pile is governed by one of the 
following considerations s 

(i) The capacity of the soil to support the load transmitted 

by the pile. 

(ii) The capacity of the pile to transmit load to the soil. 

The former consideration is generally governed by the 
properties of the adjoining soil. The later consideration 
is known as structural capacity which depends on the strength 
of pile material and the pile cross-section. The pile material 
can fail either due to excessive buckling or excessive com- 
pressive stress. The axially loaded pile foundations are 
designed for the capacity of the soil to support the load 
transmitted by the pile divided by a suitable factor of safety. 
However, if the pile is bent the load carrying capacity of 
pile is dependent on the structural strength of the pile because 
the residual stresses developed in the pile due to initial 
bending constitute significant portion of allowable stress or 
yield stress for the pile material. For economy, the straight 
piles are designed such that both the pile material and soil 
are stressed to their maximum allowable value. If the economi- 
cally designed straight piles get bent during driving 4 the load 
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that can be imposed on such a pile without overstressing the 
pile material is certainly less than the load carrying capa- 
city of the soil. In the present investigation it is assumed 
that initial bending occurs in an economically designed pile 
which causes reduction in load carrying capacity of the 
straight pile. 

The work described herein is highly idealized. This is 
deliberate because it is argued that until the general mecha- 
nics are fully appreciated and understood, there is little 
hope for improving design methods to enable the load carrying 
capacity of initially bent pile to be quantified or for improv- 
ing construction techniques to enable piles to be driven 
straight. However, the computed behaviour should be qualita- 
tively correct whereas quantitatively approximate. 

The soil is represented by Winkler medium in this inves- 
tigation. In this concept the soil pressure at a point is 
assumed to be directly proportional to the deflection at that 
point and the constant of proportionality is known as coeffi- 
cient of subgrade reaction. In other words the soil is re- 
presented by infinitesimally small closely spaced independent 
elastic springs as shown in Fig. (1.2). The works of Biot 
(1937) and Vesic (1961) deal with the validity of this assump- 
tion for analysing the problem of beams on elastic foundation. 
It can be argued that the lateral deflections in the case of 
an axially loaded initially bent piles are small and thus it 
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may be assumed that at working loads, the coefficient of 
subgrade reaction remains constant with deflection at any 
depth. However, as pointed out by Chan and Hanna (1979) 
coefficient of subgrade reaction varies with depth which has 
been accounted for in this investigation. 

In the past, as summarized in Chapter II, attempts have 
been made by Johnson (1962), Broms (1963), Kurma Rao (1975), 

Rao et al (1977) to develop simplified methods of arriving at 
allowable working or design loads assuming coefficient of 
subgrade reaction either constant with depth or linearly in- 
creasing with depth. No attempts have been made to study the 
behaviour of initially bent piles in non-homogene ous soil. 

This aspect is studied in this investigation. 

In Chapter III initially bent bearing piles are analyzed 
for two layered soil system. 

A layer coefficient is defined as - the ratio of constant 
subgrade reactions of top and bottom layer. The initial shapes 
of the bent pile for three different end conditions viz. Fix-Fix, 
Fix-Pin, Pin-Pin are represented by polynomials satisfying both 
the natural and geometric boundary conditions. Closed form 
solutions are obtained for all the end conditions considered. 

The effect of length of pile, initial pile offset, layer coeffi- 
cient, the top layer depth and the axial load on the additional 
deflections and additional bending moments are investigated. 

The load carrying capacities of the bent piles are expressed in 
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percentage load carrying capacities of the straight piles 
under similar conditions. Number of polynomials are consi- 
dered to represent the initial shape of the bent pile for 
each case and unique relations are presented between initial 
radius of curvature and the percentage load carrying capacity 
for all the cases considered. Numerical data has been obtained 
comprehensively and results are presented in nbn-dimensional 
form. 

In Chapter IV initially bent friction piles are analyzed 
by finite difference technique using Glesser*s ( 1953 ) method 
as no closed form solution could be obtained for this case. 

The effect of soil friction on the additional deflections, 
additional bending moments and load carrying capacities are 
investigated. Also effect of linear and polynomial variation 
of coefficient of subgrade reaction on the capacities of 
initially bent friction pile in non-homogeneous soil is inves- 
tigated. Down-drag effects are also investigated for few 
cases. Results are plotted in non-dimensional form. 

Finally in Chapter V, conclusions are drawn from the 
investigations and suggestions for future research on this 
topic are given. 






CHAPTER II 


LITERATURE REVIEW 


2.1 General : 

Previous investigations on tent pile may be broadly 
classified in three groups: (l) Field and laboratory observa- 
tion of pile bending, (2) attempts to assess the performance 
of bent piles and recommendations for design and (3) theore- 
tical analyses to determine the pile capacity assuming cer- 
tain initially bent shapes. By far, most of the reported 
investigation is on field observations starting about 25 years 
ago. 

2.2 Field and Laboratory Ob s ervations : 

Parsons and Wilsons (1954) reported a case where 800 
end bearing piles of composite material were driven through 
a silty stratum to bed rock. Each pile was 40 m. in length 
the lower part being a 0.27 m. diameter steel tube and the 
upper part being 0.35 m. diameter corrugated shell. The piles 
were driven in preformed holes of limited depth and a heavy 
mandrel was used for driving the first 15 m. Slope indicator 
measurements showed the greater part of deviation from the 
vertical to be in the lower section of the piles with slopes 
exceeding 15°. At least 63 piles were bent to varying degrees % 
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the bases of some piles being upto about 2 m. from their 
intended position. An approximate method of assessing the 
safe loads was put forward, but this has limited general 
application in view of the assumptions made concerning the 
lateral soil pressure acting on the pile shaft. 

According to Cummings (1956) piles are some times 
distorted into long pitch helical curves, some times into 
long sweeping bends and occassionally into sharp bends 
called dog-legs. 

Johnson (1962) made field measurements on 27 m. long 
composite piles in sand and observed that the out of plan 
position of the pile base was of the order of approximately 
10 percent of the pile length. Like Parsons and Wilsons (1954) 
he presented a method of safe load estimation that depended 
on a knowledge of the deflected shape of the pile and the load 
distribution along the pile shaft. 

Perhaps the most serious deviations from plan position 
were observed in end bearing pre-cast piles driven in Gothenburg 
through a soft clay stratum to bed rock at 60 m. depth. The bases 
of these piles deviated by upto 11 m. from their intended 
position. The deflected shape of these piles which were made up of 
several segments was a gentle curve with the slope increasing with 
depth (National Swedish Council for Building Research, 1964). 
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Worth et al (1966) observed similar trends on much 
shorter piles driven' through stiff clay and sand to a hard 
clay stratum. Also the distribution of load along the 
length of pile was measured. 

Hanna (1968) reported the results of three tests on 
long H section steel piles driven through firm clays to bed 
rock at 44 m depth. A Wilson slope indicator was used to 
determine the plan position of these piles with depth. In 
contrast to the pile bending observations of others these 
piles had two bends and the field observations are reproduced 
in Fig, (2.1). 

This diagram gives the measured components of pile 
deflections for two of these piles. In addition to these 
observations extensive load testing was carried out to arrive 
at a permissible load. Several points of practical signifi- 
cance emerged from this work. (1) From inspection of pile 
driving records it was not possible to decide if the pile 
bending had taken place. (2) Pile deviations were greater 
than normal pile to pile spacing in closely spaced groups * 

(3) When a bent pile is loaded overstressing of the pile 
material resulted. (4) The toe level of. the bent piles is 
higher than that of a perfectly straight pile and end bearing 
piles may not be at bed rock horizon, (5) The initiation of 
the pile bending is due to the generation of asymmetrical 
stresses in the pile due to eccentric pile tip reactions and 
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eccentric driving inherent to all pile driving work. The 
reverse curve is due to primarily vertical weight component 
of inclined pile forcing the pile to bend. 

There were instances in which two piles formed a con- 
junctive (U) so that one rose while that other was driven 
down and vice versa (Johnson, 1968). 

Fellinius (1971) dealt with the bending of slender 
pre-cast piles which were driven in segments of 10 to 14 m. 
in length and spliced in the field by means of rigid steel 
couplings. An inclinometer for measuring the initial shape 
of bent pile was described and case studies are given for 
number of piles. 

Fellinius (1972) dealt with the bending of slender pre- 
cast concrete piles with double curvature. 

Kim et al (1973) conducted full scale lateral load 
tests on H piles of 13.3 m length. Twenty 10BP42 steel piles 
were driven to refusal. The piles were driven in three groups 
of six piles each with two isolated single piles (one vertical 
and one batter). Each pile was instrumented with 22 strain 
gauges and a 3.81 cm^ steel tube to be used with a slope indi- 
cator. It was reported that the vertical piles deflected 
about the weak axis of the pile section during driving and 
the out of verticality was in the range of 0.15 m to 1.05 m. 

It was also observed that all piles were twisted during driving. 
The angle of twist ranged from 14° to 22° in eight piles, 26° to 
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40° in six other piles and 5° to 7° in the remaining piles. 

In general, the piles showed a tendency of twisting in the 
clock-wise direction. 

Kim et al (1974) have described a new type of inclino- 
meter, known as Digitilt-Minimose, which is a 22 mm diameter 
unit.Digitilt-Minimose is designed to be used in tubes that 
are smaller than the usual plastic or aluminium tubes, that 
accommodate the standing Digitilt Model. This new type is 
useful where multiple piles are joined by rigid caps. 

Chan and Hanna (1979) reported the results of an 
experimental study at laboratory scale on the behaviour of 
bent (dog-legged) single piles when subjected to vertical 
loads. Varying degrees of bend were adopted for both fric- 
tion and end bearing piles. The test piles which were embedded 
in cohesionless soils were instrumented to measure the distri- 
butions of axial load, bending moment and transverse shear 
along the pile shaft. Following conclusions are drawn from 
their experimental study. 

1. The load settlement behaviour of bent friction piles is 
governed by the degree of bend. 

2. A friction pile with large degree of bend has greater 
ultimate load than a straight or slightly bent pile. 

3. Repeated loading appears to have considerable influence 
on the pile behaviour, particularly with respect to residual 
loads. 
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4. Large moments and shears are locked in the pile on 
unloading that can be directly related to the non-recover able 
axial and lateral movements of the shift resulting from 
external loading. 

5 * Wien a pile is bent in one plane the axial loading of 
that pile causes two dimensional movements of the shaft. The 
tendency to undergo further bending diminishes with distance 
from the bend. 

6. The coefficient of subgrade reaction for the initially 

bent piles increases with depth but decreases with horizontal 
pile movements. 

A brief listing of above case histories is presented 
in Table 2.1. 

From the above review it may be concluded that the pile 
bending or deviations during driving does occur in practice 
and it is not restricted to a particular soil type, pile 
system, or pile length. Some of the many inter-related factors effec- 
ting this process are discussed early in Chapter I. 

2»3 Performance of Bent Piles and Design Recommendations ; 

In the second category of the reported work on bent 
piles, a few attempts have been made to provide acceptable 
rules for rejection of a pile that is out of straightness 
after driving. With hollow steel piles a rule of thumb has 
been used where by a tubular pile is rejected if part of the 
base cannot be seen when light is lowered inside it (Hanna, 1970) . 
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Norwegian Practice is to place a 50 mm diameter tube 
in the pile or welded along side of the pile in such a way 
that a small inclinometer may be lowered after driving 
(Bjerrum, 1957) and where the radius of curvature is less 
than 3000 m. the pile is usually rejected. 

Swedish Practice recommends (Heilman, 196?) no reduc- 
tion in the allowable load for concrete piles where the 
average radius of curvature of the pile axis is greater than 
400 m over a length of 10 m, or- 100 m. over a length of 2 m. 
For the steel pile the largest permissible deviation at a 
joint should not be more than 1 in 75 and the maximum allow- 
able deviation is 1/1000 of the total pile length in the 
bottom 90 percent of the pile length. The maximum allowable 
deviation in top 10 percent of the pile length is 1/300 of 
the measuring distance. 

2.4 Theoretical Analyses ; 

No significant theoretical work has been done on bent 
piles. Timoshenko (1936) was the first to obtain a solution 
for the case of a uniform pin-pin imperfect column in an 
elastic continuum which can be treated as a pile. 

Glick (1948) presented an expression for buckling load 
taking the case of an imperfect pile bent in n half sine 
waves, as sum! rig a differential equation for the elastic line 
of the pile. 



Parsons and Wilsons (1954) presented solutions for 
finding safe loads on dog-legged piles. 

The first systematic approach in the analysis of bent 
piles was presented by Johnson (1962) using Winklers hypothe- 
sis. 

Broms (1963) presented a method for calculating the 
maximum lateral deflection, the maximum soil reaction and 
maximum bending moment for an axially loaded initially bent 
pile. He has in fact considered imperfect pile. 

Marcus (as reported by Johnson, 1968) presented a solu- 
tion for pinned-pinned initially bent pile . He was the first 
to show that the residual stresses due to bending are consi- 
derable. He used Winklers hypothesis with constant sub- 
grade modulus. 

Francis et al (1965) have considered imperfections 

and loading eccentricities. They have shown that a pile 

with initial imperfections has the same buckling load as 

straight pile, however, the lateral deflections just before 

the critical load P is reached, may be large. 

cr 

Burgess (1975) has studied the directional instability 
caused by the " flutter" of a pile during driving. According 
to him the problem of a pile driven into a uniform soil is 
analogous to the aerodynamic problem of a cantilever with 
lateral elastic support, placed in a wind tunnel or flume. 
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aligned towards the flow. He predicted the critical depths 
for piles by studying the changes of the systems fundamental 
frequencies of oscillation when it is loaded by follower 
forces . 

Madhav and Kurma Rao (1975) have presented a rational 
approach for the analysis smooth uniformly bent piles sub- 
jected to axial load for different end conditions using - 
Winkler% hypothesis. They have also analysed imperfect piles 
for various end conditions. Methods for calculating the 
lateral load capacity of a bent pile have also been presented 
by them (1977) along with design charts. Rao et al (1975,1977) 
have presented solutions for cases of bent and imperfect piles. 
Kurma Rao (1975) analysed settlements of initially bent piles 
and the dog-legged piles. 

Arumugam (1977). used elastic continuum model for the 
analysis of initially bent and imperfect piles. He emphasized 
the significance of initial shape of bent pile on the load 
carrying capacity by taking number of initial shapes formed 
by the first root and second root of the transcendental 
equation for pin-pin case. He also considered the shapes 
formed by the combination of the first root and second root. 

He applied reliability analysis to the initially bent and 
imperfect piles. 

Omar (1980) presented a theoretical and experimental 
study on the directional stability of the driven piles. He 
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considered the flexural and torsional stresses and investi- 
gated the stability of the path followed by the pile. He 
pointed out that skin friction, horizontal subgrade modulus, 
the distribution of these parameters with depth and point 
resistance influence stability of the pile. 

In summary, the previous investigations on bent piles 
can be classified into three categories. In the first cate- 
gory of work Parsons and Wilsons (195^) , Johnson (1962), 

Hanna (1968), Fellinius (1971), Kim et al (1973), Chan and 
Hanna (1979) and others reported the field and laboratory 
observation of pile bending. In the second category of work 
different institutions and authorities established thumb rules 
for the rejection of already bent piles. In the third category 
of work capacity of the already bent piles were evaluated theo- 
retically by Glick (1948), Parsons and Wilsons (1954), Johnson 
(1962), Broms (1963), Kurma Rao (1975), Rao et al (1977), 
Arumugam (1977) and others. Burgess (1975), Omar (1980) and 
others have worked on the directional stability of piles. 

A review of the literature available on the initially 
bent piles suggests that no work has been done on the initially 
bent piles in non-homogeneous soils which is the situation 


often arises. 



Table'Zli Summary of reported pile bending measurements 
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CHAPTER III 


AXIAL LOAD CAPACITY OF INITIALLY BENT END 
BEARING PILES IN LAYERED SOIL 

3.1 General ; 

The effect of two layered soil system on the flexural 
behaviour and the load carrying capacity of initially bent 
pile is investigated analytically in nondimensional terms. 

A modulus of subgrade reaction is used to define the soil 
stiffness* The stiffness of the bottom layer is assumed to 
vary in some constant proportion of top layer. 

Theoretical solutions for the initially bent piles in 
soils having constant subgrade modulus along the pile length 
were given by Johnson (1962), Broms (1963) Madhav and Kurma 
Rao (1975)* In practice seldom a soil stratum is encountered 
which has a constant soil profile. To take this into account 
the present investigation includes the analysis for soil 
profile which has varying modulus of subgrade reactor-with 
depth. 

The analysis presented herein provides a qualitative 
study on the effects of soil layering on the initially bent 
end bearing piles* Two layered system is considered in this 
analysis, however, the same analysis can be extended to any 
number of layers. 
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3.2 Capacity of Straight Steel Pipe Pile ; 

A pile can fail either due to adjoining soil failure 
or due to structural failure. The capacity of pile due to 
soil, failure depends on shear strength characteristics of 
the soil and various other factor which are discussed in many 
standard Soil Mechanics books and shall not be discussed in 
this thesis. 

The problem of selecting allowable stress values: for 
piling is a very complex one. In practice two relatively 
simple approaches are used. One is to assign fixed values of 
allowable stress to pile material and other is to consider 
pile as a column. Using the concept of pile as a column is 

< 

a rational approach and most . widely adopted by practising 
engineers. In this analysis this concept is used for calcu- 
lating the allowable stresses in pile material. The design of 
a pile column is a function of six basic parameters. 

1. The yield strength of the material. 

2. The stiffness of material. 

3. The equivalent unbraced length. 

4. The amount of load P acting on the column. 

5. The eccentricity of the load. 

6. Hie cross-sectional area A of the column. 

The load carrying capacity of straight steel pile is 
governed by the allowable stress in the column. IS-800 (1962) 
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has adopted the secant formula as the basis and the factor 
of safety adopted is 1.68. A similar but simple formula is 
given by American Institute of Steel Construction. According 
to AISC specifications 1962, 


F 


[1 - (&) 2 /2C c 2 ] F v 


for — < C 

r — c 


(3.1) 


where , F, 

K 

L 

r 

C c 

E 

F. 


a 


= 10^76x10° for — > C 

(KL/r) 2 r - o 

= axial stress permitted in the absence of 

bending moment 

= effective length factor 

= actual unbraced length 

= governing radius of gyration 

= Euler ratio = |> 2 E/F y ] l/2 

= Youngs Modulus of the pile material 

KIa ( gL}3 


(3.2) 


= factor of safety =* 


5 + ± ri 

3 8C^ 


8C. 


(3.3) 


F y = specified minimum yield point for the type of 
steel being used. 


Equation (3.1) applies where the ratio Kl/r is less 
than the value of C . Equation (3.2) applies where the ratio 

W 

KL/r is greater than C c . 

For the structural steel confirming to IS 226-1962, 

F y = 2600 kg/ cm 2 , and E = 2.047 x 10^ kg/ cm 2 . In this 
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investigation the ratio of internal diameter to the external 
diameter of the pipe :pile is tentatively taken as 0*9. 

3*3 Capacity of Bent Pile : 

As with straight pile the capacity of bent pile is 
checked from soil consideration as well as structural consi- 
deration of the pile material. 

3.3*1 From Soil Considerations ; 

The capacity of bent piles in homogeneous soils depends 
on the strength characteristics of the soil material and a- 
method for estimating its capacity is given by Broms (1963). 

In this method the maximum soil reactions should not exceed 
one third of the ultimate lateral resistance. The allowable 
soil reactions are given by, 

Cohesive Soils - q a11 * 1.5 ^ (3.4) 

where q u is unconfined compressive strength. 

Cohesionless Soils- = P v ( i 1 ffn f>i ^ (3.5) 

where p is the effective over burden pressure and 0d is the 
v 

angle of internal friction of the soil as measured from drained 
triaxial or direct shear test. 

From soil consideration the allowable load on bent pile 
is given by, 



** U; w + qal1 


(3.6) 
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where P cr is the critical buckling load.. 

Up is maximum additional deflection along the length 
umax 

of the pile 

is coefficient of horizontal subgrade reaction. 
q a ]_]_ is allowable soil resistance. 

Cummings (1938) obtained critical load of a pile 
hinged at both ends. He concluded that critical load exceeds 
the crushing strength of the pile unless the soil is excessively 
soft. He gave the following formula for the critical load. 

P cr = (2m 2 + 2m + 1) (3.7) 

where m is the number of half waves of the sinusoidal curve 
into which the pile tends to buckle. 


L is the length of pile and 

El is the flexural rigidity of the pile. 

Davisson ( 1963 ) presented charts for calculating the 
buckling loads for different end conditions. 


Kurma Rao (1975) has given the following expressions 
for the critical buckling load of piles with fix— fix, fix-pin, 
pin-pin end conditions. 


For fix-fix end condition, 

2 


or 


,2 
4tc m 

z m 2 z * 2 

HI 


( 3 . 8 ) 
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For fix-pin end conditions. 


2 4z £ 

3 _ _JE . sl 

° r ‘ kz m ” 2 

m 


For pin-pin end conditions. 


(3.9) 


cr 


- (^94V 


m 


m 


(4.494)' 


( 3 . 10 ) 


where, 


m 

L 

T 

T 

El 

D 

K 

"h 


L/T 

length of pile 

[El/KD]^^ for cohesive soil. 

[El/hjJ 1 / 5 for cohesionless soil, 

flexural rigidity of the pile, 
diameter of pile, 

coefficient subgrade reaction, constant for 
cohesive soil. 

coefficient of subgrade reaction for cohesion- 
less soil. 


3.3.2 From Structural Considerations ; 

AISC does not give specifications regarding the design 
of bent piles. However, it gives the design procedure for 
structural members subjected to combined compressive and 

bending stresses as given below* 

f f f 

Wien |p^ <. 0.15 IT" * IT" — 

*a . a b 


(3,11) 
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When 


'a 


> 0.15 


C f, 
m b 


(i- -&) F b 

p ' u 

e 


< 1.0 


(3.12) 


where. 


f & = computed, axial stress 

F & = axial stress permitted in the absence of 
bending moment 

f^ = computed compressive bending stress 

F b = compressive bending stress permitted if bending 
moment alone existed 

C m = coefficient dependent on column curvature 
caused by applied moment 


10.476 x 10 6 
(Kl^/r 2 ) 


K = effective length factor in the plane of bending 
= actual unbraced length in the plane of bending 
r b = radius of gyration about the axis of concurrent 
bending 


It is assumed that the cross-sectional shape of the 
pile remains constant before and after bending. C m varies, 
from 0.85 to 1,0 and in this investigation C ffi is tentatively 
taken as 1.0. To simplify the analysis it is further assumed 
that f /F 0 is always greater than 0,15. This is justified as 
the axially loaded initially bent piles are being considered. 
Accordingly Eqn. (3*12) is reduced tos 



<Z-D 



is the flexural stress due to initial bending of the 
pile 

S 2 is the flexural stress due to additional bending of 
the pile 


is the axial stress in the pile 


3.4 Statement of the Problem: 


An idealized bent pile in a layered soil is illustrated 
in Fig* (3*1). Under the action of the axial load P there is a 
tendency for the pile to bend and deflect as shown by the dashed 
lines. The relative curvature, the magnitude of the displace- 
ments, and the number of nodes is a function of the pile length 
L, rigidity of the pile El and the soil modulus of the foundation 
material K. u^, q^, and S-^ are the lateral displacement, 
the lateral pressure, the bending moment, and the shear force 
in the pile prior to the application of load P. The pile is 
assumed to be end bearing in which the entire applied load P 
is transmitted to the tip of the pile u 2 , q 2 , and'S^ are 
respectively* the additional lateral displacement, lateral 
pressure, bending moment and shear force on the pile due to 
the addition of load P. 


The soil reaction q at any depth y is expressed as a 


linear function of the lateral deflection u at that point i.e. 
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q = k u, where k is coefficient of subgrade reaction and 
u is the displacement. This represents the assumption that 
the beam column in a soil medium may be treated analytically 
as a beam column restrained by a series of infinitely closely 
spaced independent elastic springs "with stiffness k. Vesic 
(1961) presented several examples to show the validity of 
this assumption. Fig. (3.2) shows some of the variations of 
the soil modulus with depth for different types of soils 
and the one assumed in present investigation. 

In analyzing a layered soil system the soil modulus 
of the bottom layer will be expressed in terms of the soil 
modulus of the top layer. Both relatively soft and relatively 
hard upper layers are considered. Fig. (3.3) shows the varia- 
tion of soil modulus with depth that has been assumed in the 
analysis. The soil modulus of the top and bottom layers are 
and respectively. The ratio K^/K^ is termed as layer 
coefficient R • 

3.5 Assumptions and the Governing Differential Equation ; 

The following assumptions are made in the analysis. 

1. The pile is of uniform cross-section with homogeneous 
material and is bent in a single plane with a smooth curvature 
and can be idealized as beam. 

2. The moments, shears and displacements caused by the 
forces q acting on a curved beam are determined by applying 
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the same forces to a straight beam of the same length. This 
assumption has the same effect as determining the moments, 
and so forth, in an arch as though the arch were a beam. 

The computed moments and so forth will be greater than those 
that actually exist in the structure. The assumption therefore 
gives conservative results. 

3* Failure of the pile material take place first since 
the initial stresses due to pile bending could be considerable. 
4. The axial force in the pile along the length remains 
constant. 


The first two assumption are inherent in the derivation 
of the governing differential equation for the behaviour of 
axially loaded initially bent pile as given by Johnson (1968). 
The fourth assumption is valid only for this chapter. 


Johnson (1968) and Rao et al (1975) have previously 
derived the differential equation to describe the behaviour 
of initially bent piles from equilibrium considerations for 
the cases of K equal to constant with depth and K varying 
linearly with depth respectively. Eqn. (3.14) is the governing 
differential equation -where the flexural stiffness El, and the 


axial load P are constant along the length of the pile for 


the case of constant soil modulus with depth. 

d\u d 2 u„ d 2 u 1 

El g + P g + kD u p = -P 

dy dy *■ dy 


(3.14) 
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Adopting conventional soil modulus K = kD, eqn. (3.14) can 
be written as, 

d\u d 2 Up d 2 u 1 

El £ + P f + K u„ = -P 4 (3.15) 

dy dy z z dy^ 

To make a linear differential equation K is assumed to be a 
function of y only. The resulting equation becomes a fourth 
order linear, non-homogene ous differential equation with 
constant coefficients which can be solved in the closed form. 


3.6 Closed Form Solution: 


The governing differential equation for the layered 

case can be written as, 

,4 ,2 ^2 

d u 0 . d Up, d u-, 

EI —M * P — .§£ *K t u = - P —i for 0 < y < L 

dy dy dy 


(3.16) 


EJ _2b + p —f. + u a = - P 


dV 


dy 


dy 


d y 


2 for <, y <. L 
(3.17) 


where, 


u~. = Additional deflection of the pile in the upper 

2t 

layer due to axial load P 

u„, = Additional deflection of the pile in the bottom 

2b 

layer due to axial load P 
= Soil Modulus for top layer 
= Soil Modulus for bottom layer. 
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The eqns. (3.16) and (3.17) can be written in non- 
dimensional form by making use of non-dimensional parameters 
given by Matlock and Reese (i960). 

A relative stiffness factor T is defined as 

T - r Si l 1 /^ 

“ L K t J 

Since the additional deflections are a function of pile material 
and its cross-section together with stiffness of supporting 
soil it is proposed to non-dimensionalize the pile length, 
initial deflection and additional deflection with respect to T. 

Non dimensionalized depth coefficient z = y/T 
Non dimensionalized length of pile z ffi = L/T 
Non dimensionalized thickness of top layer = z ± = L-^/T 
Non dimensionalized pile offset = C/T 
Non dimensionalized initial deflections w x = u-j/T 
Non dimensionalized additional deflections are 

W 2t = U 2t //T 5 w 2b = U 2b y/T 


Substituting above non dimensional parameters in the eqns. (3.16) 
and (3.17) they can be written in non-dimensional form as. 


+ w 2t = 


for 0 < z < z. 


(3.18) 


d 2 w 2b d 2 w-, 

l?“ + R o W 2b = "7 fOTZ l< 2 <^ 


(3.19) 
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where, 


a 



and R 

o 



Three pile boundary conditions such as fixed-fixed, 
fixed-pinned, pinned-pinned are analyzed. The initial shape, 
shear force and bending moment for each case are shown in 
Fig. (3.4). 


3.6.1 Pile with Fixed-Fixed end Conditions ; 

In dimensionless form, the initial shape of initially 
bent pile with fixed-fixed boundary conditions may be repre- 
sented by, 

w ] _ = | [1.25 + 5.0 - 8.75 + 3.5 -3 ] ( 3 . 20 ) 

2 m z m z m z m 

where C/T is the non-dimensional offset of the pile tip from the 
vertical axis. The function w^ needs to satisfy both geome- 
tric and natural boundary conditions of the problem. The geo- 
metric boundary conditions refer to the deflection and slope 
and natural boundary conditions refer to bending moment and 
shear force at the boundaries. Hoff (1956) has pointed out 
that it is generally preferable to select functions that 
satisfy both the geometric and natural boundary conditions, 
because as smaller number of terms are sufficient for an 
acceptable engineering approximation. 

Equation (3.20) satisfies explicit geometric boundary 


conditions , 
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( '"l ) Z=0 " ° > 
dw 1 

( dz~^z=z = ® 
m 


(w n ) 

1 Z=Z 


m . 


C_ 

T 


dw.. 

9 ( ~y=> ) _. = 0, 

5 ' dz z=0 5 


as well as the implicit natural boundary conditions, 

( k )„_n t 0 9 ( 


d 2 w-, d 2 w n 

"T2 K=0 t 0 » ^ ~2 \ = z ^ 0 

dz dz m 


d^w-, 


d^w-, 


V-L V/-9 U. W-j 

( "7? )z =° ^ 0 ’ ( ~3 J 


dz" z=z m 


£ 0 


Therefore, Eqn. (3.20) can be used to represent the 
initial shape of the pile with fixed-fixed end conditions. 

Let the particular solutions of Eqn. (3.18) and (3.19) 
be respectively, 


W 2t 

2 3 

= + a n z n + a 0 z + a-,z 

0 112 3 

(3*21) 

W 2b 

= b + b z + b„z 2 + b-,z^ 
o 1 2 3 

(3.22) 

W 2t 

is the additional deflection of the pile 

in the 


top layer and w^ is the additional deflection of the pile 
in the bottom layer. Substituting w 2t and w 2fc with their 
derivatives in Eqns. (3.18) and (.3*19) and comparing the 
coefficients of equal powers of z on both sides the constants 
a 0 , a 1 , a 2 , ay b Q , b^, b 2 and b^ can be obtained as. 


o 


-2.5 a C 

— * - /-\ m 


210 a £ 

~T 


m 


m 


» 5 'M, C 

Rz 2 f 
o m 


210 cc C 

c 2,'i T 
R z 
o m 
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-30 a C , 420 a 2 C 


3 T 


5 T 


, b. 


m 


m 


-50 a C 420 or C 

R z 3 T R 2 Z 3 T 
o m o m 


a. 


105 a C 

~pr t 

m 


105 a C 

R 2 ^ T 
o m 


-70 a C 
, 5 f 


70 a C 
5 T 


in 


R z 
o m 


The solution for homogeneous part of Eqns. (3*18) and. (3.19) 
are , 

+ 


w, 


2t 


sin (m-^z) cosh (n^z) + cos (m-^z) sinh (m^z) 

/ 

+ sin (m^z) sinh (m^z) + cos (m-^z) cosMnnjZ)/ 


for 0 < z < z n 


(3.23) 


w. 


2b 


= sin (m!jz) cosh (m^z) + cos (m^z) sinh (nu, z) 

+ D^' sin (m^z) sinh (m^z) + cos (m-^z) coshCm^z) 


where, 


m. 


m-, 


2 


0.5 + | 
R 


and. 


for Z-, < z < z 

1 — — m 


m 2 n a 

= 0.5 - jr- 


(3.24) 


oxl/2 ^ a . m 2 f R osl/2 a 
J and m 2 “ v-jj) - f 

of the differential equations (3*18) and (3.19) 


The total solution 
is given by # 


f If 

- w 2t + w 2t 
“ w 2b + w 2b 
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The constants C-^, C^, D^, and C^» C 2 » D^, D 2 can 
be evaluated from the following 4 boundary conditions and 
four compatibility conditions. 


At 

z = 0 

W 2t = 

0 

and dw^/dz 

= 0 

At 

z = z 1 

W 2t = 

W 2b 





dw 2t 

dw 2b 

d 2 w 2t d 2 w 2b 



dz 

dz * 

dz 2 ' 

dz 2 



d3w 2t 

„ dW 2t 

d3w 2b ^ „ 

dw 2b 




' a HlzT 

dz- 5 

dz 

At 

2 s z^ 
m 

W 2b = 

0, 

dw 2b / dz = 

0 


In the case of bent piles, three types of stresses 
namely, the residual stress , the bending stress under load 
P and axial stresses are acting and their effect is cumulative. 
The various expression for stresses are obtained by making use 
of general elastic flexure equation of beams. 

Residual stress 

S 1 = It T [ ^2 + 30 “T” 105 T* + 70 “3 3 (3.25) 

Z m Z m z m z m 

Flexural stresses in the top layer S 2 ~ | | C^ a 3 z + 2a 2 

Or O 

-sin (m^z) sinh (m 2 z) {D-^n^-m^) - 2D 2 m- ] m 2 } 

^e@S (m^z) cosh (m 2 z) {D 2 (m 2 -m^ ) + 2 D p m p m 2 } 
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-sin (m-^z) cosh (m^z) (C^m^-m-^ 2 ) - 2C 2 n^m,,} 
-cos (m-^z) sinh (m^z) {(^(m^-in^ 2 ) +• 2C-^ 

( 3 . 26 ) 

Flexural stresses in the bottom layer S 0 = § - [6b,z + 2b 0 

€L <L 1 D <d 


-sin (m^z) 
-cos (m-jz) 
-sin (m-^z) 
-cos (m-jz) 


sinh (m^z) {D-jCm^-mjj 2 ) - 2D 2 m^} 

cosh (m 2 z) {D 2 (m 2 2 -m-[ 2 ) + 2D^ 

cosh (m^z) {c^(m 2 2 -m^ 2 ) - 2C 2 

sinh (m 2 z) {C^Cm^-m^ 2 ) + 2C-[ m^}] 


(3.27) 


The expressions for S-^ and S 2 are obtained by making 
use of elastic flexure equation for beams. 


Axial stress S, = — ■ • T • 5 - (3.28) 

0 m(D -d^) 

where D is the external diameter, and 
d is the internal diameter. 


3.6*2 Pile with Fixed-Pinned end Conditions e 

In dimensionless form the initial shape for fixed- 
pinned end conditions may be represented by. 



(3.29) 

c 

where ijr is the dimensionless offset of the pile tip and 
z w is the dimensionless length of pile. 

Hi 
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Eqn. (3*29) satisfies the explicit geometric boundary 
conditions, 


r dw i 

^1^=0 = 0 » (w l>z=z = I ’ ( "3z } z=0 = 0 

m 

dw-, 


1 

as well as the implicit geometric boundary conditions ( -dth=z * 0 


d 2 w- 


m 


1 

The explicit natural boundary condition ( ) 


is satisfied. 


, ^ z=z 
dz m 


= 0 


The implicit natural boundary conditions 

2 Z 

d w n d"\yu 

0 • < —4 > z -o^ 0 ’ ( 

dz z 2_u dz^ z_u 

are also satisfied. 



z=z, 


m 


f 0 


Therefore, Eqn. (3.29) can be used to represent the 
initial shape of the bent Fix-Pin pile. 

Assuming the particular solution to be of the form, 


W 2t = 

a Q + &1 z + a 2 z 2 + a 3 z 3 

for 0 < z <_ z-^ 

(3.21) 

w 2b = 

b Q + b lZ + b 2 z 2 + b 3 z 3 

for z n < z < z„ 
1 — — m 

(3.22) 


Substituting Eqns. (3.21) and (3.22) together with (3.29) 


in the Eqns. (3.18) and (3.19) respectively and equating the 
equal powers of z on both sides the. constants a Q , a^, a 2 » a^, 
b 0 , b 1 , b 2 ,and b^ can be obtained as, 

v 2 

. «"■**» jCL'UiCL V- _ 

If 


— ot C 

„ 2 f 


282 a 4 

__ 


C . 

m 9 Oa — 


J Q 


m 


m 


R z £ 
o m 


282 a 2 0 

R 2 z 
o n 
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-45 a 

R z 3 
o m 


570 a C 

?r5“ T 

o m 


141 a G 
„ 4 T » 


d 2 = 


1KL a C 


R z 
o m 


)5^a C 


b 3= ; 95 z a 5 f 

o in 


The solution for the homogeneous part of the differential 
Eqns. (3.18) and (3.19) is given by Eqns. (3.23) and (3.24). 

The total solution of the differential Eqns. (3.18) 
and (3.19) is given by, 

w 2t = w 2t + w 2t 
w 2b = w 2b + w 2b 


The constants C-^, C^, and C^, C^, D^, can be 

evaluated from the following four boundary condition and four 


compatibility conditions. 


At z * 0 


At z = z* 


w 2t " 0 


dw 2t 

“d -# 


dw 2t dw 2b 
W 2t " w 2b» — * ~ 


d 2 w 0 , d 2 w, 
ZV _ £ 


d3w 2t . . dw 2t 


d3w 2b . - dw 2b 


At 2 = z. 


0 


0 
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Knowing the deflected shape under load P, the various 
expression for stresses can be expressed as follows. 

Residual stress, 

S 1 = |||[-^ + 45-%-l4l-4 + 95 - 3 ] (3.30) 


m 


m 


m 


m 


Flexural stress S 2 in the top and bottom layers are 
given by Eqn. ( 3 . 26 ) and (3.27) respectively. 

Axial stress . S ^ is given by Eqn. (3.28). 

3.6.3 Pile with Pin-Pin End Conditions ; 

. If a pile is just driven into a rigid stratum and its 
head is not fully encased it can be treated as a pin-pin pile. 


The initial shape for pin-pin end conditions in dimen- 
sionless form may be represented as, 

w-. = 2. [q # 5 - 7 -3 + 11 *■— r - 4,5 -3 ] (3.31) 

x x m z J z z 

mm m 

Eqn. (3.31) satisfies the explicit geometric boundary 


conditions, 

( W P 2 =0 = 0 ^ (w lU - f 

m 

and also implicit geometric boundary conditions, 

dw.. dw 1 

( “cE ) z =o ^ 0 ^ ( ”3z ) z=z m ^ 0 

The explicit natural boundary conditions. 
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d 2 ^ 

( — ) ~ o and. 

dz z z=u 

are also satisfied. 


p 

d w. 


a„2 z=z TO 
dz m 


0 


The implicit geometric boundary conditions, 

d^w-, d\r-, 

( t )„ ^ £ 0 and ( i ) £ 0 

dz" 5 z=z m 


are also satisfied. 

Therefore Eqn. (3*31) can be used to represent the^ 
initial shape of the bent pin-pin pile. 


Let the particular solution of the differential Eqns. 
(3.18) and (3*19) for pin-pin case be represented as. 


w, 


w 


2 t " “0 T " 1 " ’ " 2 " ' "3 

» . . . ? - 3 


- 2 ♦ 

3 

,2 


a Q + a-jZ + a 2 z" + a^z^ for 0 < z £ z 1 (3.21) 


2b 


* b Q + b 1 z + b 2 z + b 3 z D for < z < z ffl (3.22) 


Substituting Eqns.(3.2l) and (3.22) together with the 

initial shape (3.31) in the governing differential equations 

(3.18) and (3.19) respectively and equating the equal powers 

of z on both sides the constants a Q , a^, a 2 , a^ and b Q , b^, b 2 , 

b-, can be obtained as, 

3 


264 or C 

—TP f 9 


m 


b 0 = 


264 a 2 C 
0 2 4 T 
R o z m 


a, - 


42 a C 
T 


540 C 

~T~ f 


» . b 


m 


a 


1 


42 a C 
R 2 3 f 
R o z m 


540 a~ C 
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a. 


152 a C 

- -pr~ t s 

m 

90 a C 


m 


b 0 = 


152 a C 

~ ZT T 
R z 
o m 


90 a C 

R z 5 T 
o m 


The solutions for the homogeneous part of the diffe- 
rential Eqns* (3.18) and (3.19) are given by Eqns. (3.23) 
and (3.24). 

The total solution of the differential Eqns. (3.18) 
and (3.19) are given by. 


W 2t 

~ W 2t 

f f 

+ W 2t 

w 2b 

= W 2b 

+ w 2b 


The constants C 1 , C 2 , D ±t and C^, c' 2> d£, can be 
evaluated from the following four boundary conditions and 
four compatibility conditions. 


At z = 0 


At z =. Z- 


d 2 w. 


At z = z_ 


m 


w 2t 

= 0 

and 

w 2t 

II 

4 

w 

dw 2t 
dz ~ 

d 3 w, 

C 

It . „ dw 2t 

d3w 2b 

d2 

j + a 

dz 3 

w 2b 

« 0 , 

d 2 w 2b 

dz 2 


dz 


2t _ n 

2 0 


dw, 


2b 


p 

d w. 


tt 


d 2 w, 


2b 


dz £ 


dz' 


2 


dw. 


+ a 


2b 


dz 


= 0 
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Knowing "the deflected shape under load P, the various 
expressions for stresses can be expressed as follows? 

Residual stress 


S 1 2 T T ^ + !32 — ^ - 90 ] (3«32) 

z z 2 ^ 

m m t m 

Flexural stresses S 2 in the top and bottom layers 
are given by Eqns. (3.26) and (3.27) respectively. 

Axial stress S-^ is given by Eqn. (3.28). 

3.7 Other Possible Shanes? 


A bent pile need not necessarily assume the shapes 
as given above. In fact the pile bending and the types of 
stresses developed due to bending is a complex phenomena 
and is difficult to solve rigorously. Even for the most 
idealized case there can be number of shapes satisfying all 
the natural and geometric boundary conditions as suggested 
by Hoff ( 1936 ). Different initial shapes are assumed for the 
same length and initial pile off-set to study the effect of 
shapes on the load carrying capacity of initially bent pile. 

The shape of an initially bent pile for ary end condi- 
tions can be represented by the following generalized polyno- 
mial. 


w. 


f 


t-i 

jL Z 


+ a. 


+ a- 


m 


7 + a 4 


m 


m 


* 4 + 

— zr + 

m. 


m 


( 3 . 33 ) 
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In total six boundary conditions (two geometric 
boundary conditions and the natural boundary conditions at 
each end) enable to determine the six constants in the 
Eqn. (3*33). 

3*7.1 Fix-Fix End Conditions ; 

The geometric boundary conditions to be satisfied for 


this case are* 







dw n 


at z = 0 

Q 

It 

H 

£ 

and 

X 

— 

= 0 

at z = z 

m 


and 

dw^ 

= 0 

The 

natural boundary conditions 

to be 

satisfied 


d^w-, 


d\/^ 


at z * 0 

■ — rt £ 0 

dz d 

and 

dz 5 

f 0 


2 

d Wt 


d^w-, 


at z = z m 
m 

dz^ 

and 

1? 

t 0 


Assuming arbitrarily the shear force at top and bottom 
of the bent pile as 6F ] _ and 6F 2 , (for mathematical convenience) 
and substituting the geometric boundary conditions in Bgiii (3.33) » 
six simultaneous equations can be formed which can be solved 
to obtain six constants as given by, 

10-1. 5F 1 + 0.5F 2 

a o } a-^ = o ; a 2 = 4 “ 

-10-3* 5F-j -1. 5F ? F l +F 2‘ l4 

^ f a 5 = 4 
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For each set of F-^ and F^ values such that 20-3F-^+F2 £ 0 
and 20+F^-3F2 £ 0 an initial shape satisfying the natural and 
geometric boundary conditions for fix-fix end conditions can 
be formed; 


3 * 7 * 2 . Fix— Fin End Condition; 


The geometric boundary conditions to be satisfied ares 


at z = 0 

H* 

If 

o 

and 

dw^ 

at z = z 

m 


and 

dw-^ 

“3F 


The natural boundary conditions to be satisfied are: 
at z = 0 


d 2 w^ d^w-, 

£ 0 and 4 £ 0 


at Z at z 


'm 


dz 

d 2 w-, 


dz* 


= 0 


and 


dz - 

d^w-j 

dz' 


/ 0 


Assuming the shear force at depth K equal to zero and 
the slope at z = equal to © and substituting these 
conditions together with geometric boundary conditions in 
Eqn. (3.33). Six simultaneous equations can be formed which 
can be solved to obtain six constants as given by, 

(-40+6) (30K 2 -15)-6Q 


a 0 ~ ® » 


a l = 0; 


60K Z -72K z + 18 


^5 

a 3 


20 - (12K Z - 8) (-4© + 6) 

60K 2 - 72K, + 18 

4b 4b 

( 0-2 ) — * -3a^ f 3^2 




43 


For each set of K z and Q values such that a 2 0, 
a 3 £ 0, a^ + 4a^ + 10a^ f 0 an initial shape satisfying the 
natural and geometric boundary conditions can be formed. 

3.7.3 Pin-Pin End Conditions ? 

The geometric boundary conditions to be satisfied ares 


dw. 


at z = 0 

o 

it 

H 

> 

and 

VJL VV ~|^ 

Hz 

/ o 

at z = z„ 
m 

W = £ 

I T 

and 

dw-^ 

" dz 

yL 0 


The natural boundary conditions to be satisfied ares 



d 2 w-j_ 



d^w. 


at z = 0 

= 0 

and 

1 

dz^ 

jL 0 


dz 2 




d 2 w-. 



d^w. 


at z = z_ 
m 

i. 

= 0 

and 

i 

t o 

dz 2 



' dz^ 


Assuming the slopes at the top end and bottom end of 
the piles 0-^ and 0^, respectively and substituting the geo- 
metric boundary conditions together with the boundary condi- 
tions of zero bending moment at the ends in Ean. (3.33) results 
in the formation of six simultaneous equations which ban be 
solved to obtain the six unknowns as given below. 
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3.8 Initial Curvature and Capacity of Bent Piles ; 

Hanna (1968) highlighted the effect of initial radius 
of curvature on the initial bending stresses in the bent pile. 
Kurma Rao (1975) has shown the effect of R/D ratio on the 
percentage load carrying capacity of bent pile using Winkler 
model where R is the initial radius of curvature and D is the 
diameter of the pile, Arumugam (1977) found using the elastic 
continuum model, the relation between R/D ratio and the 
percentage load carrying capacity of bent piles for pin -pin 
end conditions for different configurations to be unique. 

Different configurations are formed by the combination 
of the first two roots of the transcendental equation. In this 
investigation, the effect of R/T ratio on the percentage load 
carrying capacity of bent steel pipe piles having different 
configurations is studied for fix-fix, fix-pin and pin-pin end 
conditions using Winkler model. 


Prom calculus the expression for radius of curvature 
can be written in non-dimensional form as, 


R 

T 



The value of minimum radius of curvature along the length 
of the pile is chosen and the corresponding value of the capacity 
of bent steel pipe pile is evaluated as explained previously. 



3.9 Results and Discussion; 


To study the effect of axial load on additional deflec- 
tions (wg) and bending moments (Mg) the parameters z^ = 5.0, 

R q = 1*0, ;=r = 0.01 and ^ =22.0 are kept constant and the load 
on the pile is varied, a is the load carrying capacity of 
straight pile under similar conditions. Figs. 3* 5a, 3.5b and 
3.5c show the variation of additional deflections (wg) for fix- 
fix, fix— pin, pin-pin end conditions respectively. Figs. 3.6a, 
3.6b and 3.6c show the variation of additional bending moments 
(Mg) for fix-fix, fix-pin and pin-pin end conditions respec- 
tively, It is observed that for all the end conditions the 
deflections and bending moments are directly proportional to 
the load on the pile. For fix-fix end conditions the maximum 
additional deflection (wg) occurs at a depth of 1,25 and 3.75 

simultaneously where as the maximum additional bending moments 
occur at the ends and also at depths of 1.25 and 3.75. For 
fix-pin end conditions maximum addition deflection (wg) occurs 
at a depth of 3*75 whereas maximum additional bending moment 
(Mg) occurs at the top end and at a depth of 3.75 simultaneously. 
For pin-pin end conditions both maximum additional deflection 
and bending moment occur at a depth of 3.75.' 

To study the effect of initial off-set on the additional 
deflections (wg) and bending moments (Mg) the parameters z ffi = 5.0 
R *» 1.0, a = 0.75a r , and T/r = 22,0 are kept constant and 

v O 

C/T is varied. Fig. 3.7a, 3.7b, and 3.7c show the variation of 
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additional deflections (wg) for fix-fix, fix-pin and pin-pin 
end conditions respectively. The corresponding figures in 
Figs. 3.8a, 3.8b and 3.8c show the variation of additional 
bending moments (Mg) for the respective end conditions. As 
shown in all the figures the additional deflections (wg) and 
additional bending moments (Mg) are directly proportional to 
initial off-set of the pile. The positions of maximum deflec- 
tions and maximum bending moments are same as for previous 
case. 

To study the effect of layering on the variation of 
additional deflections (wg) and bending moments (Mg) the 
parameters z m = 5.0, z ± = 2.5, C/T = 0.01, and T/r = 22.0 are. 
kept constant and R Q is varied. Figs. 3.9a, 3.9b and 3.9c show 
the variation of additional deflections (wg) with the change 
in strength of bottom layer (R q ) for fix-fix, fix— pin and 
pin-pin end conditions respectively. Figs. 3. 10a, 3.10b, 3.10c 
show the variation of additional bending moments (Mg) with R Q 
for fix-fix, fix-pin and pin-pin end conditions respectively. 

As expected when the bottom layer is weaker than the top layer 
(R = 0.25) the bending moments (Mg) and additional deflections 
(w 2 ) increases. The increase in additional deflections (wg) 
and bending moments (Mg) is higher in the bottom weaker layer 
compared to top stronger layer. In contrary when the bottom 
layer is stronger than the top layer 2(R Q = 4.0) the deflections 
(Wg) and bending moments (Mg) decrease all along the length of 
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pile and the maximum reduction is found to be in stronger 
layer. 

The capacity of few bent piles from the soil failure 
considerations are calculated using the formula given by 
Broms (1963) and buckling loads given by Kurma Rao (1975) . 

A typical example is shown below. 

For a steel pile of length 11.22 m, external diameter 
30 cm, internal diameter 27 cm with an initial off-set of 2.24 cm 
embedded in a homogeneous cohesive soil with an undrained shear 

O 

strength of 0.34 kg/cnr the nondimensional parameters can be 
obtained as follows. 

Moment of inertia I = ^ (30^ - 27^) = 13673.73 cm^ 

Area of cross-section A = -7“- (30 2 - 27 2 ) = 134.3 cm 2 

Radius of gyration r = (l/A)“^ 2 = 10.09 cm 

Terzaghi (1955) has given the relation between coeffi- 
cient of horizontal subgrade reaction . .. (t/ft^) and undrained 

O 

shear strength (t/ft ) as K = 33 Cu. 

This can be written as, 

K = 1.1 Cu 

where, K is coefficient of horizontal subgrade reaction in 



2 

Cu is undrained shear strength in kg/cm . 
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. * . Relative stiffness factor 


/.vil _ /. 2.047 x10 6 x13673.73 

4V KD ~ 4 Cf.lxO'.34) 30 " 


223.5 cm 


The non-dimensional parameters can be obtained as, 

z m = 5.02, -jjj- = 0.01 , = 22.15 

~ 5.00 ~ 22.00 

The details of calculations for obtaining the capacity of 
bent pile from soil failure considerations are shown in Table 
3.1. All the parameters shown are dimensionless. 


It can be seen from the above example that the load • 
carrying capacity of bent pile from soil failure criteria is 
much higher than the capacity calculated from structural 
failure criteria. Hence in this investigation only structural 
capacities of the bent piles have been obtained. 

Fig. 3.11a, 3.Hb and 3.11c show the variation of load 
carrying capacity with z ffi for fix-fix, fix-pin and pin-pin 
end conditions respectively. The parameters R Q = 1.0 and 
T/r = 22.0 are kept constant. It can be seen from the figure 

# 

that for C/T = 0*01 there is no significant change in the pile 
capacity after z ffi = 5.0. However, for higher off-set values 
longer piles provide higher load capacities. 

Fig* 3.12a, 3* 12b and 3.12c show the relation between 
load, c arry ing capacity of bent pile with C/T for fix— fix, fix- 
pin and pin«*pin end conditions respectively. The parameters 
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R 0 = 1*0 and T/r = 22.0 are kept constant. It is observed 
in all the figures that the load carrying capacity reduces 
very rapidly for short piles with increasing off-set. 

Fig. 3.13a, 3.13b and 3.13c show the relation between 
C/L, initial off-set expressed in terms of pile length, against 
capacity of bent pile for fix-fix, fix-pin, pin-pin end condi- 
tions respectively. At smaller pile offsets the behaviour is 
not well defined but at higher offsets the capacities of 
longer piles are larger than shorter piles. 

Fig. 3.14a, 3.14b and 3.14c show the effect of slender- 
ness parameter T/r on the percentage load carrying capacities 
of bent piles with fix-fix, fix-pin and pin-pin end conditions 
respectively. The parameters z m = 5.0, R Q = 1.0 are kept 
constant. At larger offsets percentage load carrying capacity 
increases as the slenderness parameter increases. This increase 
is attributed to lower load carrying capacities of straight end 
bearing piles. However, for smaller offsets, there appears to 
be a slenderness parameter at which capacity will be maximum. 

Fig. 3.15a, 3.15b, 3.15c show the effect of R Q on the 
load carrying capacity of bent pile with fix-fix, fix-pin, 
and pin— pin end conditions respectively. The parameters 

z m “ 5 *°» Z 1 “ 2 * 5 °» C//T = °* 01 ^ T//r = 22 ’ 0 are kept COns ~ 
tant. As seen in these figures for all end conditions the 

load carry ing capacity increases with the increase of R 0 due to 

increase in the strength of the bottom layer. However the 
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magnitude of increase is ansi g inf icant. This may be due to 
high Young's Modulus of steel piles. If the initial pile off- 
set is high the initial beading stresses would be high and 
hence allowable load and consequently the additional deflec- 
tions would be small. If the Initial offset is small the 
axially load capacity of b eat pile would certainly be high 
but the eccentricity will b e smaller and thus additional 
deflections and bending stresses would be small. Thus in both 
cases the small additional deflections would make the soil 
strength an insensitive parameter to the load carrying capa- 
city of bent pile. For this reason we find that the effect of 
the strength of soil is insignificant on the load carrying 
capacity of bent steel piles. 

Figs. 3.l5d, 3.15e, 3.15f show the effect of the depth 
of top layer ‘ on the lead carrying capacity of bent pile for 
fix-fix, fix-pin and pin-pin end conditions respectively. The 
parameters ^ . 5.0, C/1 . 0.01 and T/r = 22.0 are kept eons- 
tant. The parameter is given values of 0.25 and 4.00 and 
z 1 is varied between 1.25 and 3.75. It can be seen from these 
figures that for all end conditions when R q = 0.25 (i.e* when 
the bottom layer is weaker) idncrease in the thickness of top 
layer increases the load carrying capacity where as when 
R 0 * 4*6 increase in the thickness of top layer decreases the 
load carrying capacity. However, the change In load carrying 
capacity with change in thiefeness of top layer is insignificant. 
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The effect of initial offset on the additional deflec- 
tions of bent piles at the allowable load on the bent pile 
is shown in Fig. 3.16a, 3.16b, 3.16c for fix-fix, fix-pin, 
and pin-pin end conditions respectively. In this figure the 
additional deflections are expressed as a percentage of initial 

deflections for various values of z„. It is seen that at 

m 

larger offsets the percentage change in additional deflection 
is smaller than at smaller off-sets. For example, in fix-pin 
case for z ffl = 5.0 a change in the value of C/T from 0.02 to 
0.04 (i.e. 50 percent increase) reduced the percentage addi- 
tional deflections from 1.0 percent to 0.4 percent (i.e. 60 per- 
cent reduction) in additional deflections). This clearly shows 
that the effect of initial bending stresses becomes predominant 
as C/T increases. For pin-pin end conditions the maximum 
additional deflections are negative. However, the increase in 
the magnitude additional deflections reduces at higher offset 
values. Fig. 3.17a, 3.17b, 3.17c show the additional bending 
moments as a percentage of initial bending moments for the 
same pile with fix-fix, fix-pin and pin-pin end conditions 
respectively. 

It is recognized that the load carrying capacity of 
bent pile depends on the initial shape of the pile. To study 
this aspect number of shapes have been considered and analyses 
were conducted keeping initial offset constant, for fix— fix, 
fix-pin, pin-pin end conditions. For all shape| .pid 
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conditions considered the parameters z =5.0, R =1.0, 

T/r = 22.0 are kept constant and C/T is varied between 0.01 
to 0.05. Three typical shapes, shape I, shape II, shape III 
for each end condition are shown in Fig. (3.18). The load 
carrying capacities for shape I, shape II and shape III are 
evaluated with different initial offset values for all t h e 
end conditions and are shown in Fig. (3.19). 

It is interesting to note that for all shapes and 
different offset values, a unique relationship is found between 
R/T and load carrying capacity for each end condition (Fig. 3. 20). 
For all end conditions the load carrying capacity increases as 
radius of curvature R/T increases. At a particular radius of 
curvature the capacity is maximum for fix— fix end conditions 
and minimum for pin-pin end conditions. Similar curves between 
R/D (Ratio of radius of curvature to diameter of the pile) and 
load carrying capacity were presented by Kurma Rao (1975), 
Arumugam (1977) which is independent of pile end conditions. 

In their analysis it is assumed that the allowable stress in 
the pile material remains constant irrespective of end condi- 
tions and length of the pile. However, in the present inves- 
tigation the steel pipe pile design is based on the column 
concept in which allowable stresses changes with pile length 
and end conditions. For this reason different relations are 
obtained for different end conditions. 



Table 3.1: Capacity of Bent Pile from Soil Failure Consideration. 
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load applied TO PILE 

DIAMETER OF pile 

additional soil pressure 
mobilised due to p 
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CHAPTER IV 


INITIALLY BENT FRICTION PILES IN LAYERED AND 
NON-HOMOGENEOUS SOILS 


4*1 General: 

Pile foundations derive their strength both from end 
bearing resistance and skin friction. In Chapter III, the 
analysis of end bearing piles has been presented. In this 
chapter the influence of friction along the pile on the 
load carrying capacity of the bent pile is studied by modify- 
ing the governing differential equation for both, cohesive 
and cohesionless soils* Three sub-soil conditions as shown 

in Fig* (4.1) are analyzed, (l) Two layer system of cohesive 

non-horn ogeneous 

soils, (2) Two layered system of cohesionless soil, C 3 V soils 
with polynomial variation of subgrade reaction. The polyno- 
mial variation of subgrade reaction takes into account any 
deviation of soil modulus from the idealized variations such 
as soil modulus constant with depth and linearly varying with 
depth (Matlock and Reese, I960). The effect of down drag is 
also studied. Because a closed form solution is difficult 
to obtain for the above mentioned models finite difference 
technique is used based on the method suggested by Glesser 

( 1953 ). 
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In the case of cohesive soils, the coefficient of sub- 
grade reaction is assumed to be invariant with depth for 
each layer (Terzaghi,1955) and thus the total load taken by 
the skin friction is directly proportional to the surface area. 
This will provide a linear distribution of axial load in the 
pile. In the case of cohesionless soils the coefficient of 
subgrade reaction is assumed to be linearly varying with depth 
(Terzaghi, 1955) and since the skin friction Is a linear func- 
tion of depth, the axial force distribution becomes parabolic. 
In the case of polynomial variation of coefficient of subgrade 
reaction also the axial load distribution is assumed to be 
parabolic. 

4.2 Finite Difference Technique : 

4.2.1 Basic Concept ; 

In this method the continuous bent pile Is divided into 
N segments. The unknown function, for the present case addi- 
tional deflection, is replaced by N algebraic variables at 
the discrete points.. Consequently, the differential equations 
which describe the behaviour of continuous initially bent piles 
are replaced by N simultaneous algebraic equations in these . 
variables. In the process the derivative of the unknown func- 
tion at the point is approximated by an algebraic expression 
consisting of the value of the function at that point and 
several near-by points. Hence, the closer the points are to 
one another the better is the agreement between the derivative 



56 


and its algebraic approximation and the more accurate will 
be the solution. However, as the number of points increases 
so does the number of simultaneous equations that must be 

solved. 

\ 

4.2.£ Solution Procedure? 


The pile is divided into n segments. In addition to 
this two segments at the top of the pile and two segments at 
the bottom are imagined as shown in Fig. 4,2. The following 
relationships can be written about node i in difference form. 


dw 0 

(-*£). = 
v dz '1 


( w ?> i+1 - 
_ 2 h 


,2 r 

d w 


dz' 


2 >, - 


d\r 


dz 


2 ), 


5 


- 2(w ? ) i 1 (v, 2 ) i-i 
' h 2 


(w 2 ) 1+2 - 4 ( w 2h+l * 6 ^ w 2h ~ 4 t vf 2 h-l' l '^ W 2^ 


i<*2 


h 




(4.1) 


where , 


h = increment length 
(w 2 ) i « the additional deflection at node i. 

Substitution of the above differentials in the governing 
differential equation results an algebraic equation with the 

(w 2 ) 1+2 , (w 2 ) ltl , (w 2 ) ± , (* 2 >u ^zh-2 as 

This type of equation can be obtained for all nodes (0 to n) 

which forms n + 1 equations in n + 5 unknowns. 
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Therefore, four boundary conditions are needed to get the 
complete solution. Two boundary conditions at top and two 
boundary conditions at bottom for fix-fix, fix-pin, pin-pin 
are provided as discussed in previous chapter. 

A convenient method of solving these large number of 
simultaneous equations "has been suggested by Glesser (1953) • 
The steps in Glessers method are as follows. 


1. A recursive equation is set-up for deflections at each 
point along the pile, starting at the top. 

2. The coefficients in the recursive equation are computed 
along the length of the pile. 

3. The deflection at the bottom of the pile is calculated 
from recursive equation and the boundary conditions at the bottom 
of the piles. 

4. The deflections for all upper points, starting at the 
point next to the bottom and progressing upward, can then be 
computed using the recursive equation. 



The governing differential equation for a friction 


bent pile is given, by, 

El + §7 E p ( y ) §7)3 + K 7 (y) = 0 (4,2) 

dy y 

where, u = u^ + u 2 and u- L and u 2 are initial and additional 
deflections respectively . 
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The solution of Eqn. (4.2) depends upon the form of 
variation of axial force and also on the variation of 

Soil modulus The most simple case for which closed 

form solution is available is that when the axial force as 
well as the soil modulus are invariant with depth. The analysis 
of the bent piles with friction is presented for. the cases of 
cohesive and cohesionless soils. 


Since the finite difference method is applicable only 
for continuous functions, the discontinuity in the soil modulus 
at the interface of two layers is approximated as shown in 
Fig. 4,3. 

4.4 Analysis of Bent Friction Piles in Cohesive Soil s 


Assuming the soil modulus to be invariant with depth 
in each layer and axial lop/.' in pile to be linearly decreasing 

with depth Eqn. (4.2) reduces to. 


El 


El 


.4 
d u 


2t 


3£A 1^2t P ?1 du 2t . 

T / 9 ~ T Hrr + 


"dy* + P (1 - ?1 L' dy 2- L dy 

d 2 u 


L 

?1 

4- P — 

2 r L 


K t U 2t 


du. 


dy 

for 0 < y <, 

d\i~ L n y-L 1 

— f + p(1 -'’ir-' , 2-r> 

dy 

L l 

+ K b u ?b - — 


_1 

dy 


(4.3) 


dy 


2b p _2 
2 ” r L 


0 du 


2b 


dy 


y- L i 


9c 


d 2 u-j 

dy 2 


cp 0 du, 

+ P ¥ "dy for L 1 - y - L 



59 


where <p ^ and 9 2 are "the skin friction factor in the top 
layer and the bottom layer respectively. All other terms are 
explained in the previous chapter. 

The relative stiffness factor T is defined as T = [fl] 1 ^ 4 

Non-dimensionalising deflections and length by using T as 

Non— dimensionalized depth coefficients z = y/T 

Non-dim ensionalized length of pile z m = L/T 

L 1 

Non-dimensionalized thickness of upper layer z-^ = 


n 

Non-dimensionalized pile offset = 
Non-dimensionalized initial deflections w-^ = 



2 

Non-dimensionalized additional deflections w^ = -7p 
Non~dimensionalized total deflections w = 


Substituting the above non-dimensional parameters in 
the Eqns. (4.3) and (4.4) they can be written in non-dimensional 
form as, 


d 4 w 0 . d 2 w p dw , 

— — ^ + a (1-p.jZ) — ”2 - “Px dz + w 2 t 

dz dz 


d 2 w. 


dw. 


= ^(l-Pq 2 ) T ^1 “dz 


— ** -j 

i + ocp, for 0 < z < z.. (4.5) 


,4 
d w, 


fp + a(l-p 1 z 1 -p 2 (z-z 1 )) f y “ a P 2 “dz + R o w 2b 


d 2 w. 


2b 


dw. 


2b 


dz 


d 2 w n 


dw. 


U. I»n JL 

-a(l-{3-^Z2~ P2^ z “’ z l^ ' , 2 + a ^2 " Sz 

\Jl£* 


for z ± < z < z m 


(4.6) 
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where, 


and R 


The remaining terms are explained previously. 

For convenience the total pile length is divided into 

n elements such that a node falls at the interface of two 

layers. The node n-, at the interface of two layers is given 
z 

n i = -r x n 

m 

z 

The length of pile increment is given by h = . 


Since the initial shape of the bent pile w^ is an 
assumed polynomial it can he differentiated and substituted 
in Eqn. (4.5) and (4.6). At any general node i the right hand 
side of the differential equations (4.5) and (4.6) can be 
written as, 

2 

f d w 1 dw-, 

RHS (I) = -ad-f^zd)) (q4>z=zCl) + “ P 1 ( "c#z=z(I) 


for 0 < i < n n 


(4.7) 


RHS'(I) 


V*. VV-j 

(1-P 1 z 1 + P 2 z x - p 9 z(l)) («— ~) 


£^2^ z=z ( i) 


+ a P 2 <-a?>z«(i) £or n i i 1 i 


(4.8) 


where z(l) is the depth of the node i which can be obtained 
as, z(l) = (i) (h). 
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Substituting the difference equations (4.1) in Eqns. 

(4.5) and (4.6) the left hand side of these equations at any 
nod .0 1 can foa wniirfcGn as y 

[2W2(l+2) + FA(I) W2(I+1) + FB(I) ¥2(l) 

+ FC(l) W2(I-1) + 2W2(l-2)] (4.9) 

where W2(l) is the additional deflection, w 2 at node i, 

FA(I), FB(I), FC(I) are the coefficients of ¥2(1+1), W2(l), 
and W2(l-l) respectively which can be obtained as 
below: 

FA(l) = -8.0 + 2ah 2 (l-p^zd)) - a^h 3 

FB(I) = 12.0 - 4ah 2 (l-p^d)) + 2h 4 

FC(l) = -8.0 + 2ah 2 (1-^ z(l)) + ap^h 3 

for 0 <_ i <_ n-^-5 

FA(I) = -8.0 + 2ah 2 (l-P-jZd)) - ap-jh 3 

R —1 

FB(l) = 12.0 — 4ah 2 (l— P-^z(l)) + 2h 4 (l+ — (i+5— n-^)) 

FC(l) * -8.0 + 2ah 2 (l-p^zd)) + ccp.^ 3 

for xij-5 <. i ■<_ 

FA(I) = -8.0 + 2ah 2 (l-P 1 z- L + p^ - p 2 z(l>) - ap^ 3 

FB(I) » 12.0 - 4ah 2 (l-p^ + p^ - P 2 z(l)) 

R —1 

+ 2h 4 (1+ (i+5-r^)) 

FC(l) •« -8.0 + 2ah 2 (1-p-^ + p^ - P 2 z(l>) + ap^ 3 

for n^ < i < n-^ +5 
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* 


FA(I) = -8.0 + 2ah 2 (l-p + P 2 z i " p2 2 ^ 1 ^'” a ^2 h ^ 

FB(I) = 12.0 - 4och 2 (l-P^ + § 2 z ± - P 2 z(l)) + 2h 4 R c 

FC(I) = -8.0 + 2ah 2 (l-p^ + - p 2 z(l)) + ap ^ 3 

for n^+5 <. i <. n 

Defining, 

RHS (I) = 2h 4 RHS ’(I) (4.10) 

"the differential equations (4.5) and ( 4 , 6 ) at any node i can 
be written in terms' of the additional deflections of neighbour- 
ing nodes as, 


2W2 (1+2) + FA(I) ¥2(1+1)+ FB(I) W2(l) + FC(I) ¥2(1-1) 

+ 2W2 (1-2) = RHS (I) (4.11) 

4,4.1 Fix-Fix End Conditions ; 

The initial shape satisfying both natural and geometric 

boundary conditions for fix-fix end conditions is given by, 

r 2 3 4 5 

w 1 = f [1.25 ~ + 5.0 - 8.75 -Tj- + 3.5 -3 ] (3.20) 

Z m z m z m z m 

Differentiating Eqns. (3.20) and substituting in Eqn. 

(4.7) and (4.8) gives, 


RHS '(I) 
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RHS (I) = -a(l-§ lZ + p z - (3 z(l)) | [ ^4 

z ^ 

m 

+ 30 241 _ 105 + 70 2.(1 1- ] 

2 m 

+ ap 2 2. [2,5 dil + io d.J ) 2 _ 35 S i. ? 3 » 

2 Z 2 

m m m 

+ 17.5 ] for n 1 <_ i <_ n (4.13) 

z J j_ — 

m 

RHS (I) = 2h 4 RHS'(I) 

Boundary conditions at top: 

At z = 0, w 2 = 0 and ~~&£ ~ ® 

Writing these boundary conditions in finite difference form. 


W2(l) - W2(-l) = 0 

2 h 

gives, W2(l) = W2(-l) 

Writing the governing differential equation in finite difference 
form Eq. (4.11) at node 1, and substituting top boundary condi- 
tions W2(0) = 0 and ¥2(1) = W2(-l) leads to, 

- 2.0 + A Fb|i) W2 ^ + 2.0 + FB(l) ¥2 ^ + 2.0+FB(l) 

* 

= A(l) W2(2) + B(l) W2(3) + C(l) 

Similarly writing the governing differential equations in 
finite difference form Eq. (4.11) at node 2 and substituting 
W2(0), ¥2(1) leads to 

¥2(2) = A(2) W2(3) + B(2) ¥2(4) + C(2) 
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where , 

A(2) = 

B(2) = 

C(2) = 


r fa( 2) 

+ B(l) 

fc(2)! 

FB(2) 

+ a(i) 

FC(2) 


-2.0 


FB(2) 

+ A(l) 

FC(2) 

RHS(2) 

- QCi) 

FC(2) 

FB(2) 

+ A(l) 

FC(2) 


For any general node i, the recurrence relation can be 
written as, 

W2(I) = A(I) ¥2(1+1) + B(I) W2(I+2) + C(I) (4.14) 

where , 


A(I) 


FA(I) + B(I-1) FC(l) + 2A(l-2) B(I-l) I 

FB(I) + ACI-1) FC(l) + 2A(I-l) A(l-2) + 2B(I-2; 


(4.15) 



C(I) 


RHS(l) - C(I-1) FC(l) - 2A(I-2) C(l-l) -2C(I-2 
FB(I) + All-l) FC(lJ+ 2A(I-l)A(I-2)+ 2B(l-2 

(4.17) 


Writing the boundary conditions at the bottom end of the pile 
in finite difference form. 



leads to W2{n+1) = W2(n-l). 
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Writing the general recurrence relationship Eqn. (4.14) 
at node (n-l) and substituting W 2 (n) = 0 and W 2 (n+l) = W 2 (n-l) 
the deflection at node (n-l) can be obtained as, 

W2 < n -« ■ rifei] 

Knowing W2(n-1), W2(n) and the recurrence coefficients 
A(l), B(l), C(l) at all nodes, by the process of back substi- 
tution, additional deflections W2(l) at all nodes can be 
evaluated. 


4.4.2 Fix-Pin End Conditions: 


The initial shape satisfying natural and geometric 
boundary conditions corresponding to fix-pin end conditions 
is given by. 


w. 


• x = § [0.5 ~^2 + 7 .5 - 11.75 + 4*75 - 3 - ] 


m 


m 


z z 

m m 


(3.29) 

Differentiating Eqn. (3.29) and substituting in Eqn. (4.7) 
and (4.8) gives, 

RHS'(I) = -aCMSjzd)) § [ + 45 ^ - 141 


m 


m 


+ 95 ^ 


] + ap 1 | [ + 22.5 


m 

,2 


m 


m 


m 


- 47 

z 


2 (?)- + 23.75 z ^ I g ” ] for 0 < i £ n 


m 


m 


1 

(4.18) 
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RHS' (I) 


RHS (I) 

Writing Eqn.(4.11) at node 1 and substituting the top 
boundary conditions expressed in finite difference form, 

A(l), B(l) and C(l) can be obtained. Similarly writing 
Eqn. (4.11) at node, 2 and substituting W2(0) and W2(l) leads 
to the evaluation of A(2), B(2) and C(2). The remaining 
recurrence coefficients A(l), B(I), C(l) for all nodes can 
be obtained using Eq. (4.15) , (4.16) and (4.17). 

Expressing the bottom boundary conditions, W2(n) = 0 

d 2 

and (_ — £.) . = 0, in finite difference form and substituting 

dz^ ~^m 

in the general recurrence relation Eqn. (4.14) at node n-1 the 
deflection at node (n-1) can be obtained as, 

Knowing W2(n), W2(n-l) and all other recurrence coeffi- 
cients, by the process of back substitution, the deflections at 


= -a (l-p 1 z 1 + $2 z i_ ~ P 2 z(l)) I t ~^2 


m 


+ 45 _ i4i + 95 £(I Y -j 

V 


+ ap 2 | [ + 22.5 _ 47 


m 


m 


m 


+ 23.75 ] 


for 


n-^ £ i <_ n (4.19) 


m 


= 2h 4 RHS ' (I) 
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other nodes can be evaluated. 

4.4.3 Pin-Pin End Conditions : 

The initial shape satisfying natural and geometric 
boundary conditions corresponding to pinned-pinned end condi- 


tions is given by, 


f [1.5 z . 


sXil. _ 7,0 z ^\ 


■ini. 


+ 11.0 


4.5 


m 


m 


m 


z(I ) 5 ] 
— 5 - J 
. z m 

(3.31) 


Differentiating Eqn. (3-31) and substituting in Eqn. (4.7) 
and (4.8) gives, 

2 

RHS 1 (I) = -a (1— p-^zCl) ) [-42 + 132 ^ i 


m 


m 


RHS 1 (I) 


- 90 — ^ [ 2 ™ - 21 ^ 3 " 


m 


m 


m 


-\4 


+ 44 - 22.5 — ^ ] for 0 < i < n x 


m 


m 


(4.20) 


-a(l-P 1 z 1 + P 2 Z 1 ~ P 2 z(l)) | f" 42 ^3 

m 


+ ,132 


2(1 y 


z(i) ] + ap ^ [ 


90 


C r 1*5 


2 T L z. 


m 


m 


m 


_ 21 + 44 - 22.5 


m 


m m 

for n-, < i £ » 


(4.21) 
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Writing the top boundary conditions in finite difference form* 

W2(0) = 0, ( d2w 2) _ Mil - 2W2(0 ) + W2,( r l) _ 0 

^77 ; z=0 " 2h . 

z 

gives W2(l) = -W2(-l). 

Writing Eqn. (4.11) at node 1, and substituting 
W2(0) = 0 and W2(l) = -W2(-l) the coefficients A(l), B(l) 
and C(l), can be obtained as, 

- FA(1) _ -2.0 

-2.0 + FB(TT » ^ } ~ -2.0 + EB(l) * 

RHS(l) 

-2.0 + FB(1) 

Eqn. (4.11) at node 2 and substituting W2(0), 
coefficients A(2), B(2) and C(2) can be obtained as, 

[FA(2) + B(l) FC ( 2 ) ] 

FB(2) + A(l) FC(2) ’ 

- 2.0 

FB(2) + A(l) FC(2) 

RHS(2) - C(l) FC(2) 

FB(2) + A(l) FC(2) 

The recurrence coefficients at any general node i are 
given by Eqn. (4.15)» (4.16) and (4.17). 

Writing the bottom boundary conditions in finite 


A(l) = 

C(l) = 

Writing 
and W2(l) the 

A(2) = 

B(2) = 

C(2) = 


difference form, 
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Cl 

W2(n) = 0 , (—-£-) = W2(n-fl) -2¥2(n) + W2(n-l) _ Q 

V, 2 z=z ,2 “ 

d z m h 

gives, W2 (n+1) = -W2 (n-1) 

Writing the general recurrence relation (4.14) at node 
(n-1) and substituting bottom boundary conditions the deflec- 
tion at (n-l)th node can be obtained as, 

W2(n - 1 > - 

Knowing W2 (n), W2(n-1) and all the recurrence coeffi- 
cients, by the process of back substitution the deflections 
at all the nodes can be evaluated. 


4.5 Analysis of Bent Friction Piles in Cohesionless Soil ; 

The assumption of soil modulus constant with depth 
in each layer may not be necessarily true in all cases. For 
normally consolidated clays and especially for granular 
soils the assumption of soil modulus varying linearly with 
depth is more appropriate. When the soil modulus varies 
linearly it is reasonable to assume the variation of axial 
load in the pile to be parabolic. For this assumption 
Eqn. (4.2) reduces to, 


d 4 u 


2 da . 2P? j. du 2i 

■r— + y -i p Y 1 2 

^ X j ^ rhr L 


EX + p(l-<p ^j) 

dy 4 X L dy 


dy 


“hi? 


u 


2t 


2 ddi- 

-P (1- 9 i ^2 ) 
jL 


_JL 2P9i X. 

dy 2 l 2 ay 


for 0 < y < 


i 


(4.22) 



70 


d 4 u 


El 


2b 


dy 




Pd_ 9 b! - » y 2 * y i 2 ~ 2yy i • d2u 


1 T 2 


-9 


) — » 
dy 


y n du„ 

- 2P<p, - -f ) -g^ + y u a 


TO ' 0 

d L ■ if 


= -P(l- 9 


'i 9 2 yi 


1 L 2 


L 


2 ^2 


y 2 -2yy x d 2 U]L 


) 


dy 4 


v Y - ] du, 

+ 2P 92 ( X- - ) -X for L x < y < L 


(4.23) 


where 9 ^ and 9 ^ are skin friction factors for cohesionless soils 
for the top layer and bottom layer respectively. 

■^ht’ ^b are coe ffioients of modulus variation for 
the top and bottom layer respectively. 


The remaining terms are explained previously. 

The relative stiffness factor T is defined as T = [ XL. ] 1 / 5 

n ht 

Non-dim ensionali zing the parameters with T as shown in 


previous sections, the following relations can be written, 


u. 


m 


— mTo ■ Z-. — 


9 » — 


_i 

T ’ 


_ ^2 

Wp y "W 


Substituting the above non-dimensional parameters in the 
Eqns. (4.22) and (4.23) they can be written in non-dimensional 
form as, 


d 4 w. 


2t 


dw. 


2 

0 d w 0 . ^■"'0+ 

l * « a-Pi z ) —¥ - 2a ih z — - w st z 

dz^ dz 2 

« d w-, dw, 

= -a (1-P-jZ 2 ) f + 2a i^l z “dF for 0 < z < z ± 

± dz 


(4.24) 


t-3|P 
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d^w 


2b • « 2 „ _ 2 „ , 2 _ .. d %b 


dz 


TT + aCl-p^ - (3 z -T - Mz 2 -2zz-.)) 


2 1 ~ r 2' 

dw 2b 


' 1 ” 'a?' 


+ (2 h z i - 2 h z) -af - •*• R 0 ’"ab 2 




-a(l-f 1 z 1 - P 2 Z 1^ “ ^ z “ 2zz l^ “JT 

dz 


dw^ 


+ 2ap 2 (z- Zl ) -gi for Z;L < z < z m (4.25) 


where , 


P 


<Pl 9 o 

1 ~ p &nd P p = o 9 

z_ z 


m m 

The remaining terms are explained previously. 


R = S* 

0 ”ht 


At any general node i the right hand side of the above 
differential equation can be written as, 

« o d2w i 

RHS (I) = -a (l-p^z (ir) ( "^■) 2=Z ( I ) 

dw 1 

+ 2ap 1 z(l) ( ar >z-z(i) lor °^ 2 i (4 - 26) 

2 

O 0 O d W 1 

RHS (I) = -a(l-p 1 z 1 ^-p 2 z 1 " ; - p 2 (z z -2z(l) z-J) z=z(l) 

dW-, 

+ 20 C§ 2 (z(l) - Zj) ( ‘"^2 ^ z=z (i) 

for z-j < z <, Zjjj (4.27) 

Substituting the difference equations (4.1) in Eqns. (4.24) 
and (4.25) the left hand side of the equation at any node i can 

be written as, 



(4.28) 


where s 


[2¥2 (1+2) + FA(I) W2(l+1) + FB(I) W2(l) 
+ FC(I) W2(I-1) + 2W2 (1-2)] 


FA(I) 
FB(I) 
PC (I) 

FA(I) 

FB(I) 


-8.0 + 2ah 2 (l-p lZ (l) 2 ) - 2ap i h 3 z(l) 
12.0 - 4ah 2 (l-p iZ (l) 2 ) + 2h 4 z(l) 
-8,0 + 2ah 2 (l-p lZ (l) 2 ) + 2ap 1 h 3 z(l) 


for 0 < i < n^-5 

-8.0 + 2ah 2 (l-P lZ (l) 2 ) - 2ap i h 3 z(l) 
12.0 - 4or,h 2 (l-P lZ (l) 2 ) + 2h 4 (Xn.j-5) 


R o 2(^+5) - z(n.j-5) 

-- I5 - 


(i+5-n-^)] 


FC(l) 


FA(I) 


FB(I) 


FC(l) 


-8.0 + 2ah 2 (l-§ lZ (l) 2 ) + 2ap 1 h 3 z(l) 
for n^-5 < i £ n-^ 

-8.0 + 2ah 2 (1-P^ 2 - p^ 2 - p 2 (z(l) 2 -2z(l) Zl )) 
- 2ap 2 h 3 (z(l) - z-^) 


12.0 - 4ah 2 (l-p 1 z 1 2 - P 2 Z 1 2 ~ ^(zd^-Szd^)) 

R z(ru+5) - z(n-,-5) 

+ 2h 4 [ z(n r 5)+ — - — — — — (i+5-r^) 


-8.0 + 2ah 2 (l-^i 2 ! 2 - $2 Z 1* “ P 2 (zd)^-2z(I) Zl )) 



for i n-^+ 5 


+ 
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FA Cl) = -8.0 + 2ah 2 (l-p^ 2 - p^ 2 _ P 2 ( z (l) 2 -2z(l) Zl ) ) 
- 2ap 2 h 3 (z(l)- Zl ) 

FB(I) = 12.0 - 4ah 2 (l-p^ 2 - p^ 2 _ P 2 ( z (l) 2 -2z(l) 2l ) ) 
+ 2h 4 R o z(l) 

FC(l) = -8.0 + 2ah 2 (l-p^ 2 - p 2 2 - p 2 ( z (l) 2 -2z(l)z 1 ) ) 
+ 2ap 2 h 3 (z(l)- Zl ) 

for n-^+5 < i < n 

Defining RHS (I) = 2h 4 RHS f (I) 

the governing differential equations (4.24) and (4.25) can be 
written at any node i in terms of additional deflections of 
neighbouring nodes as, 

2W2(l+2) + FA(I) W2(I+1) + FB(l) W2.(l) + FC(I) W2(l-l) 

+ 2W2 ( 1-2 ) = RHS (I) (4.29) 

The differential Eqns. (4.24) and (4.25) can be 
solved to find the additional deflections at all nodes for 
fix-fix, fix-pin and pin-pin end conditions exactly on similar 
lines of cohesive soil. The only difference between these two 
cases is the values of RHS (I), FA(I), FB(I), FC(l) differ 


one another 
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4 . 6 Analysis of Bent Friction Piles in Non-Homogeneous Soils ; 

In this section the effect of uniformly varying soil 
with the polynomial variation of soil modulus and parabolic 
variation of axial load in the pile is studied. The governing 
differential equation for this case takes the form: 



(4.30) 


The relative stiffness factor T is defined as. 


m r Si 1 !/(n+4) 

i - L Kl J 

Non-dimensionalizing the parameters with T as shown in previous 
sections, the following relations can be written. 

Z 


T U 1 U 2 U 

z = ¥ > z m = T’ W 1 = T ’ w 2 = " • w 


T * 


Substituting the above non-dimensional parameters in Eqn. (4.30) 
it can be written in non-dimensional form as. 


d H w 


r , d 2 W 0 “«p , n . 

£ ♦ a(l-pz 2 ) —4 -2af3z -gf + ir Q +z ) w 


dw^ 


dz 


% 


dz p 

dw n dw-, 

= -a (1-pz 2 ) — £ + 2apz (4.31) 

dz 


where « 


O . T ... 

jj ss w 


and 


r o = 


K q T 

El 


m 


The remaining terms are explained previously. 
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4.6.1 Solution s 

At any general node i, the right hand side of the 
differential equation (4.31) can be written as, 

i o d 2 w-, 

RHS (I) = -a Cl-MD 2 ) ( -a^) M(I) 

dw-, 

+ W 2 ™ ' "3z >z=z(I) 

Substituting difference equations (4.1) in equation (4.31) the 
left hand side of the differential equation about node i in 
finite difference form can be expressed as, 

[2W2 (1+2) + 'FA(I) W2(l+1) + FB(I) W2(l) 

2h 

+ FC(I) ¥2(1-1) + 2W2 (1-2)] 

where , 

FA(I) = -8.0 + 2(xh 2 (l-(3z(l) 2 ) - 2apz(l) h 3 

FB(I) = 12.0,:- 4ah 2 (l-pz(l) 2 ) + 2h 4 (r Q + z n ) 

FC(I) = -8.0 + 2ah 2 (l-pz(l) 2 ) + 2aph 3 z(l) 

Defining RHS (i) = 2h 4 RHS '(I), 

the governing differential equation at any node i can be 
expressed in terms of additional deflections of neighbouring 
points as, 

2W2(I+2) + FA(I) ¥2(1+1) + FB(l) ¥2(1) + FC(l) W2(l-l) 

+ 2W2(l-2) = RHS (I) 

The above equation can be solved for Fix-Fix, Fix-Pin, and 


Pin-Pin end conditions exactly on similar lines of cohesive 
soils with friction. 
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^*7 Load Carrying Capacity of Bent Friction Piles ; 

Finite difference technique used to solve the differen- 
tial equations for various cases mentioned in the previous 
sections provide the additional deflections under axial load- 
ing. The procedure to obtain the load carrying capacity from 
additional deflections is as follows. 


The calculated additional deflections are added to the 
initial deflections to get the total deflections, w = w^ + w^. 
From the finite difference expression for bending moment at 
node i. 


,2 

d w 


w (i+l) ~ 2w (i) 


i/,r 


n-D 


dz 


h £ 


where, w ^ represents the total deflection at node i 
the bending moments at all nodes are calculated. 


The criteria for design is. 


C f, 
m b 


(1 + F a , ) F fe 
e 


< 1 


(3.17) 


The various terms in the above equations are explained in 
Sec. 3.2. For friction piles in layered cohesive soil, 

^ P 


-j- (l-P-jjz) for 0 < z < 


~ - P 2 ( z “ z l^ for z i i 2 i z , 


m 


where , 


<Pl 


m 


9 2 


m 


and 
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For friction piles in layered, cohesionless soil, 
fp. = (1-3 -,z^) for 0 <, z <. z-j_ 

= ~f~ (1 “Pi 2 } 2 “ ^2 for z l 1 2 1 z m 

where, 



For friction piles In non— homogeneous soil, 

f a - Hr < 1 - P 22 ) 

where, p = — 

2 

m 

The highest axial load which causes the axial and 
bending stresses along the length of bent pile satisfying 
equation (3. 17) is determined by iteration process and Is 
expressed as percentage capacity of straight end bearing pile. 

^•*8 Analysis of Bent Piles Subjected to Down Drag s 

When a pile is driven through compressible strata to 
a relatively firm material, there is a tendency for the com- 
pressible strata to consolidate by an amount in excess of the 
displacement of the pile under applied load. In such a case 
the soil around the pile tries to hang up on the pile, exert- 
ing a downward drag. The drag may be of considerable magnitude 
and consequently the load transferred to the bottom portion of 
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"bent pile is increased. The analysis of bent piles with 
the drag force on the entire length of the pile was studied 
by Kurma Rao (1975) and Arumugam (1977). In this investi- 
gation the drag force factor is varied in both top layer 
and bottom layer. The solution is same as that of friction 
piles with the exception that the skin friction factor 
acquires negative .sign. 

4.9 Results and Discussion ; 

To check the validity of the finite-difference method 
a comparison is made between closed form solution and finite- 
difference solution with respect to load carrying capacity. 
Table 4.1 shows the load carrying capacities of a smooth 
bent pile with fix-fix end conditions in a homogeneous soil 
for different initial pile offsets. The parameters z m = 5.0, 
R q = 1.0, T/r = 22.0 are kept constant and C/T is varied 
between 0.01 and 0.05. Table 4.2 shows the load carrying 
capacities of a smooth bent fix— fix pile in a layered soil. 
The parameters z m = 5.0, z-^ = 2.5, C/T = 0.01 and T/r = 22.0 
are kept constant and R q is varied between 0.25 and 4.0. As 
seen from these tables the finite-difference method gives 
consistently slightly higher values, however, this difference 
is insignificant as the maximum error is 6 percent in the 
worst case. 

To study the effect of friction on the variation of 
additional deflections and bending moments the parameters 
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z m = 5.0, R q = 1.0, C/T = 0.01 and T/r = 22.0 are kept 
constant. The friction factor is varied between -1.0 to 
+1.0. Fix— fix end conditions are considered for the pile. 

In cohesive soils, the effect of friction on the variation 
of additional deflections and additional bending moments 
is shown in Fig. 4.4a and 4.4b respectively. In cohesion- 
less soils the effect of friction on the variation of addi- 
tional deflections and additional bending moments is shown 
In Fig. 4.5a and 4.5b respectively. For both types of soil 
it is seen that in general as the friction increases the 
additional deflections and the additional bending moments 
decreases. But the amount of reduction in deflection for 
cohesionless soils Is about 50 percent of that in cohesive 
soils. In both cases the reduction is more in the bottom 
part of the pile than in the upper one. Effect of negative 
skin friction (9 = -l) is more severe in both cases. 

To study the effect of friction on the percentage 
load carrying capacity of an initially bent pile having 
different lengths In cohesive and cohesionless soils, the 
parameters C/T = 0.01, R Q = 1.0, T/r = 22.0 are kept cons- 
tant, z m is varied between 2 to 6 and the friction factor 
is varied between —1.0 to 1.0. Figs. 4.6a and 4.6b show the 
variation of percentage load carrying capacity with for 
fix-fix and fix-pin end conditions respectively for cohesive 
soils. Figs. 4.7a and 4.7b show similar graphs for 
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cohesionless soil. It is seen from these figures that as 
the friction increases the percentage load carrying capacity 
increases and for negative skin friction the load carrying 
capacity decreases. In all the cases, for the same friction 
factor, the increase in the load carrying capacity is higher 
for longer piles. The effect of negative skin friction is 
more severe in fix-fix pile than in fix-pin pile. For pin- 
pin pile with the shape under consideration the load carry- 
ing capacity is 100 percent of the corresponding end bear- 
ing straight piles i.e. there is no change in load carrying 
capacity due to Initial bending. This is because the maximum 
bending stress occur at a depth of 70 to 80 percent of pile 
length and at that depth the addition of bending and axial 
stresses are less than the allowable stress for a load on 
the pile head corresponding to 100 percent capacity of end 
bearing pile. 

To study the effect of friction on the percentage load 
carrying capacity of an initially bent pile having different 
initial pile offsets in cohesive and cohesionless soils the 
parameters = 5.0, R q = 1.0, T/r =22.0 are kept constant, 

C/T is varied between 0.01 to 0,05 and friction factor is 
varied between -1.0 and +1.0. Figs. 4.8a and 4.8b show the 
effect of friction on the percentage load carrying capacity of 
fix— fix and fix-pin piles respectively for cohesive soils. 

Figs, 4.9a and 4.9b show similar graphs for cohesionless soils. 
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It is observed from these figures that friction significantly 
increases the load carrying capacity at smaller pile offsets. 
As the offset value increases the effect of friction decreases 
and approximately at C/T = 0.04, it becomes negligibly 

small. 

To study the effect of friction on the capacity of 
bent piles in cohesive soils at different T/r values, the 
parameters z m = 5.0, R q = 1.0, C/T = 0.01 are kept constant 
and T/r is varied between 10 and 26 . Figs. 4.10a and 4.10b 
show the effect of slenderness parameter on the capacity of 
bent piles for different friction factors for fix-fix and 
fix-pin end conditions respectively. It can be seen from 
these figures that the load carrying capacity attained 
certain maximum value between T/r = 15 and 20 and then gra- 
dually reduces. This is because, for highly slender mem- 
bers, for a given offset even though the initial bending 
stresses are small the additional bending stresses due to 
imposed load are large making the capacity of bent piles 
small. For bulky sections the initial bending stresses for 
the given offset are high thus reducing the capacity of bent 
piles. For cohesionless soil since practically there is no 
increase in the capacity with friction It is not shown. On 
the same reasons as mentioned above pin-pin end conditions 
are not considered. 
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In order to study the effect of layering (R ). on the 
capacity of bent friction piles the friction factor for % 
the top layer is kept constant (9 ^ 158 0 * 5 ) and the friction 
factor for the bottom layer ( 9 2 ) is varied as shown in 
Table 4.3. 

Figs. 4.11a and 4.11b show the effect of R Q on the load 
carrying capacities of fix-fix and fix-pin piles respectively. 
The parameters z ffi = 5.0, z-^ = 2.5, C/T = 0.01 and T/r = 22.0 
are kept constant and R Q is varied from 0.25 to 4.0. In 
both figures, it is seen that as R Q increases i.e. when the 
strength of bottom layer increases the load carrying capacity 
also increases. However, the percentage increase is very 
small (less than 0.2 percent) and for all practical purposes 
it can be neglected. This may be explained from the fact 
that the increase in load carrying capacity due to friction 
is compensated by the reduced resistance of the soil because 
of the reduced additional deflection. Figs. 4.12a and 4.12b 
show similar graphs for cohesionless soil for fix-fix and 
fix-pin end conditions respectively. 

Figs. 4.11c and 4. lid show the effect of thickness of 
top layer (z-^) on the load carrying capacities of fix-fix 
and fix-pin piles respectively when R Q = 0.25 (i.e. when 
the bottom layer is weaker than the top layer) and R Q = 4.0 
(i.e. when the bottom layer is stronger than the top layer). 
Similar graphs are drawn in Figs. 4,12c and 4.12d for I 
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cohesionless soil. It is seen from all these figures that 
when the bottom layer is weaker than top layer increase 
in depth of top layer increases the capacity of bent pile 
and vice-versa. 

Figs. 4. lie and 4. Ilf show the effect of negative 
skin friction of the upper layer for fix-fix and fix-pin 
end conditions respectively. The parameters . = 5.0, z-^ = 2.5> 
R q = 4.0, C/T = 0.01, T/r = 22.0 and - 0.5 are kept cons- 
tant and 9 ^ is varied. Similar graphs are drawn for para- 
bolic variation of negative skin friction in Figs. 4.12e and 
4.12f. In all the cases, It is seen that as the negative 
skin friction factor increases the load carrying capacity 
reduces. While there Is a reduction of 5 percent in the 
load carrying capacity in case of soils having linear varia- 
tion of negative skin friction there is only 1 percent 
reduction in the case of soils having parabolic variation of 
negative skin friction. 

The effect of polynomial variation of subgrade modulus 
on bent friction piles Is shown in Fig. 4.13. For fix-fix and 
fix-pin end conditions friction factor Is taken as 0.5 and 
for pin-pin end conditions friction factor is taken as 0.0 
for the reasons mentioned previously. Figs. 4.13a , 4.13b 

and 4.13c show the variation of capacity of bent pile for 
r Q = 0.25 and varying n for fix-fix, fix-pin, pin-pin end 
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conditions respectively. Figs. 4.13d, 4.13e and 4.13f 

show the variation of capacity for n = 0.75 and for varying 

r values. The effect of increase in either n or r is to 
u • o 

increase the capacity of bent pile. However, for practical 
purposes the increase can be neglected. 
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Table 4.1s Comparison of closed form and finite 

difference solution in homogeneous soil. 


% 


c 

T 

Capacity in percentage 

Percentage 

Closed Form 

— GA x 

Finite Difference 

0.01 

71.57 

71.75 

0.28 

0.02 

52.30 

52.61 

0.59 

0,03 

37.04 

37.50 

1.24 

0.04 

23.75 

24.34 

2.48 

0.05 

12.08 

12.82 

6.12 

Table 4.2s 

Comparison of 

closed form and finite 


difference solution in layered cohesive soil. 


C apacity in p e rcentage Percentage 


"o 

Closed form 

Finite Difference error 

0.25 

71.52 

71.69 

0.24 

0.50 

71.54 

71.71 

0.24 

0.75 

71.56 

71.73 

0.24 

1.00 

71.57 

71.75 

0.25 

2.00 

71.61 

71.79 

0.25 

4.00 

71.64 

71.82 

0.25 

Table 4.3s 

Variation of cp 

2 wlth V 
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4.13 EFFECT OF POLYNOMIAL VARIATION OF SOIL MODULU 

ON CAPACITIES 




CHAPTER V 


SUMMARY AND CONCLUSIONS 


It is desirable to drive the piles straight, however, 
it is seldom possible in practice to drive a pile straight 
without deviations from the vertical axis. Many investiga- 
tors have contributed to the analysis of deformed piles in 
homogeneous soils which are briefly reviewed in Chapter II. 

In this investigation an attempt is made to study 
the effect of layered soils on the load carrying capacity 
of bent steel pipe pile taking into consideration important 
governing parameters such as boundary conditions, initial 
shape, offset, soil modulus etc. 

Polynomial expressions satisfying natural and geome- 
tric boundary conditions are selected arbitrarily for fix-fix, 
fix-pin and pin-pin end conditions. The soil is idealized 
by the Winkler medium and the layer coefficient is used to 
define the ratio of the coefficient of subgrade reactions of 
each layer. 

The initial stresses in the pile are evaluated from 
the assumed initial shape using general elastic- flexure 
equation. The beam column equation is used to obtain the 
axial load response of the bent pile. 
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For the case of two layered soil system having cons— 
tant soil modulus in each layer a closed form solution is 
obtained using two boundary conditions at each end of the 
pile and four compatibility conditions at the interface of 
the two layers. The maximum bending stress in the pile is 
obtained by adding the bending stress due to additional 
deflection caused by the imposed load to the initial bending 
stresses developed due to deformed shape of the pile. The 
load carrying capacity of straight pile is determined from 
the column concept as recommended by AISC Code of Practice. 

The load carrying capacity of the bent pile is obtained 
using the interaction formula given by AISC for columns 
subjected to both compressive and bending stresses. The 
highest axial load which causes the axial and bending stresses 
along the length of bent pile satisfying interaction formula 
is determined by iteration process and is expressed as a 
percentage of the capacity of straight end bearing pile. 

All the linear parameters are non-dimensionalized by a 
relative stiffness factor which accounts the rigidity of the 
pile and the soil modulus. In the column concept of pile 
design the capacity of straight pile depends on the slender- 
ness ratio of the pile. To study the effect of slenderness 
ratio on the capacity of bent pile, a slenderness parameter 
is defined as the ratio of relative stiffness of the pile 
to the radius of gyration of the pile. 
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The load carrying capacities of bent piles with 
different initial shapes for each set of end conditions are 
determined and the initial shape of the pile is found to be 
most important parameter which determines the capacity of 
bent pile* 

The beam column equation is suitably modified to 
obtain the axial load response of bent piles with side 
friction in cohesive and cohesionless layered soils. In 
cohesive soils the coefficient of subgrade reaction is 
assumed to be constant with depth and the variation of axial 
load in the pile is assumed to be linear. In cohesionless 
soils the coefficient of subgrade reaction is assumed to be 
linearly varying with depth and the variation of axial load 
in the pile is assumed to be parabolic. To account for any 
deviations from soil having constant and linear variation 
of soil modulus with depth, the effect of polynomial varia- 
tion of soil modulus with parabolic variation of axial load 
in the pile is also studied, on the load carrying capacity 
of initially bent piles. Solutions to the initially bent 
friction piles are obtained using finite difference tech- 
nique. 

The following conclusions can be drawn from the para- 
metric study conducted on initially bent steel pipe pile 
in this investigation. 
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End Bearing Piles ; 

1. For the initially bent piles the load carrying capa- 
cities calculated from soil failure criteria are much higher 
than the capacities calculated from structural failure cri- 
teria. Thus the bent pile design is governed by the struc- 
tural considerations. 

2. In the layered soil additional deflections and addi- 
tional bending moments due to axial load decreases as the 
strength of the bottom layer increases. The load carrying 
capacity of a bent pile increases as the strength of bottom 
layer increases. In case the bottom layer is weaker than 
top layer the load carrying capacity of bent pile increases 
with the increase in the thickness of top layer. Wien the 
bottom layer is stronger than the top layer load carrying 
capacity decreases with the increase in the thickness of 
top layer. 

3. As the load on the initially bent pile increases the 
additional deflections and additional bending moments increases. 
In homogeneous soils, for fix-fix case, the maximum additional 
deflections occur simultaneously at a depth of 1/ 4th and 

3/' 4th pile length. For fix-pin and pin-pin end conditions 
the maximum additional deflections occur at a depth of 
3/4th pile length. For fix-fix pile the maximum additional 
bending moments occur simultaneously at the ends of the pile 
and at depths of 1/4 and 3/4 pile length. For fix-pin pile 
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maximum additional bending moment occurs simultaneously 
at top end of the pile and at a depth of 3/4th pile length. 
For pin-pin pile. the maximum additional bending moments 
occur at a depth of 3/4th pile length. 

4. As the initial offset increases the initial bending 
moments, additional deflections and additional bending 
moments due to applied load increases for all the end condi- 
tions. The load carrying capacity of bent pile reduces with 
the increase in the initial pile offset for all the end 
conditions. For the smaller offsets the maximum percentage 
load carrying capacity is obtained with slenderness para- 
meters between 13 and 20. At other slenderness parameters 
the capacity will be smaller. However, at larger offsets 
the percentage load carrying capacity of highly slender 
piles is large. 

5. At constant offset the load carrying capacity of a 
bent pile increases with the increase in length of the pile 
upto a certain length after which there is no signifi-cant 
increase in the capacity with the increase in length of the 
pile. 

6. The load carrying capacity of bent piles is primarily 
a function of initial shape. Irrespective of initial pile 
offset the percentage load carrying capacity of bent pile 
increases as the radius of curvature of the initial shape 
increases. For the same radius of curvature the percentage 
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load carrying capacity is maximum for fix-fix end condi- 
tions and minimum for pin-pin end conditions. 

Friction Piles ; 

1. The soil friction significantly reduces the additional 
deflections and additional bending moments. For the same 
tip resistance the reduction in additional deflections and 
bending moments due to applied load is more in cohesive 
soils than in cohesionless soils. The additional deflections 
and bending moments are increased due to negative skin fric- 
tion. The percentage load carrying capacity of a bent friction, 
pile is more than that of an end bearing pile. The load carry- 
ing capacity of a bent pile with negative skin friction is 
less than that of an end bearing pile. 

2. The effect of friction is more significant for longer 
piles than for shorter ones. 

3. The effect of friction is more significant at- small®*’ 
pile offsets than at larger pile offsets 

4. The load carrying opacity of a bent friction pile is 
insensitive to ±k« strength of bottom layer. Similarly the 
polynomial variation of soil modulus has insignificant effect 
on the "Load carrying capacity of bent pile. 

Suggestions for Further Research ? 

From this study several new ideas have emerged and 
a few of them are listed below which can be persued further. 
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1. Analysis of bent piles taking torsional stresses 
into account. 

2. Analysis of bent piles subjected to axial load, 
lateral load and moment. 

3. Analysis of pile groups with some piles bent. 

4. Drivability of bent piles. 
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